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ABSTRACT 


For  this  thesis,  an  atmospheric  propagation  code  named  ANCHOR  (Atmospheric  NFS 
Code  for  High  Energy  Laser  Optical  pRopagation)  was  developed  and  utilized  to  study 
the  propagation  of  high  energy  lasers  in  various  atmospheric  conditions  and  for  numerous 
laser  configurations.  The  ANCHOR  code  accesses  existing  industry  databases  to  obtain 
relevant  optical  properties  for  various  atmospheres  and  then  uses  scaling  laws  to  simulate 
laser  propagation  through  the  defined  environments. 

ANCHOR  accounts  for  the  effects  of  atmospheric  diffraction,  turbulence, 
platform  jitter  and  thermal  blooming  on  the  laser  beam,  and  outputs  on-target  irradiance 
and  power-in-the -bucket  profiles  for  a  wide  range  of  laser  wavelengths.  Several  known 
physical  trends  associated  with  laser  propagation  will  be  reproduced,  and  the  results  will 
be  compared  to  the  industry  accepted  propagation  code  Wavetrain. 

The  results  of  ANCHOR  studies  will  indicate  that  the  100  kW-class  high  energy 
laser  can  effectively  engage  slow-moving  targets  at  ranges  greater  than  five  kilometers  in 
clear  weather  by  delivering  enough  energy  to  melt  0.1  liters  of  one  millimeter- thick 
aluminum  aircraft  skin  in  five  seconds.  For  hazy,  turbulent,  and  rainy  conditions,  the 
laser  can  effectively  engage  targets  from  ranges  closer  than  three  kilometers,  but 
reasonable  dwell  times  are  only  achieved  for  ranges  closer  than  two  kilometers. 
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I.  INTRODUCTION 


The  United  States  Navy  is  heavily  invested  in  the  researeh  and  development  of 
high-energy  lasers  (HELs)  onboard  various  surface  and  air  platforms  for  offensive  and 
defensive  capabilities  [1],  [2],  Such  weapons,  while  extremely  complex  in  design  and 
expensive  to  manufacture,  hold  many  advantages  to  currently-employed  conventional 
weapons  systems.  From  speed-of-light  payload  delivery,  to  an  extremely  deep  magazine, 
to  engagement  flexibility  for  soft  and  hard  kills,  laser  weapons  are  often  touted  as  a 
significant  part  of  the  future  of  naval  warfighting  [3]. 

Unlike  conventional  weapons,  however,  laser  weapons  are  uniquely  affected  by 
weather  conditions  of  the  environment  in  which  they  are  employed  [4],  Effective 
development  of  these  weapons  is  closely  tied  to  the  proper  modeling  and  understanding 
of  the  role  that  atmospheric  properties  play  in  hindering  or  aiding  laser  propagation. 
Unfortunately,  simply  relying  on  experimental  studies  in  order  to  fully  parameterize  laser 
behavior  in  various  weather  conditions  is  neither  time-efficient  nor  cost-effective  [4]. 
Physics-based  computer  modeling  is  therefore  essential  to  guide  HEE-based  research. 

To  this  end,  a  streamlined  and  robust  simulation  code  has  been  developed  by  the 
NFS  Directed  Energy  Group.  The  code  utilizes  existing  atmospheric  databases  and 
governing  scaling  equations  to  accurately  and  efficiently  model  laser  propagation  in 
various  environments  for  numerous  laser  configurations.  The  code  has  been  used  to 
explore  atmospheres  with  variable  turbulence  profiles,  visibilities  and  precipitation,  as 
well  as  laser  effectiveness  with  various  beam  diameters,  output  powers,  wavelengths,  etc. 

The  first  part  of  this  thesis  will  discuss  the  history,  technology  and  advantages  of 
directed  energy  weapons.  Next,  the  thesis  will  explore  the  physical  principles  of 
atmospheric  propagation,  followed  by  atmospheric  modeling,  and  then  the  simulation 
methods  utilized  to  carry  out  the  laser  propagation  studies.  Easily,  the  thesis  results  will 
be  discussed  and  analyzed,  with  conclusions  tying  together  the  quantitative  and 
qualitative  aspects  of  the  research. 
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II.  DIRECTED  ENERGY  WEAPON  OVERVIEW 


A,  DIRECTED  ENERGY  INTRODUCTION 

A  directed  energy  (DE)  weapon  emits  energy  directed  at  a  target  in  order  to  inflict 
damage  without  the  use  of  a  projectile  [5].  The  energy  harnessed  by  a  DE  weapon  is 
usually  in  the  form  of  electromagnetic  radiation,  which  includes  microwave  weapons  and 
lasers,  but  also  can  include  particle-beams  and  sonic  weapons.  Currently,  DE  weapons  of 
military  interest  utilize  a  laser  (light  amplification  by  stimulated  emission  of  radiation)  to 
generate  and  deliver  a  concentrated  spot  of  light  at  a  specific  wavelength  and  of  high 
enough  irradiance  on  the  surface  of  a  target  to  cause  damage.  The  wavelength  range  of 
the  laser  light  typically  lies  in  the  visible  to  infrared  regime  (~0.5  pm  to  20  pm),  and  the 
energy  propagates  as  a  narrow  beam  at  the  speed  of  light  to  reach  the  target.  Damage  is 
primarily  inflicted  by  the  high-energy  laser  “by  light  dwelling  on  the  target  for  an 
extended  period  of  time  to  accumulate  deposited  energy  and  inflict  thermal  damage”  [4]. 

Easer  lethality,  defined  as  “the  capability  of  a  weapon  system  to  render  a  target 
non-functional,”  [4]  greatly  depends  on  the  target  material  the  amount  of  power 
transmitted  to  the  target.  Thermal  damage  is  considered  to  be  the  main  lethality 
mechanism  for  lasers,  resulting  in  the  burning,  melting,  penetration  or  structural  integrity 
failure  of  the  target  [4].  Eor  typical  targets,  such  as  the  carbon-carbon  ceramic  nosecone 
of  a  surface-to-air  missile  (SAM),  an  estimated  12.1  kilowatts  (kW)  of  power  over  a 
dwell  time  of  4  seconds  is  required  to  ablate  a  depth  of  1.0  cm  on  the  target  [6].  While  a 
laser  emitting  such  an  output  power  level  may  not  appear  to  be  powerful  enough  to 
neutralize  a  missile  traveling  at  Mach  3.0,  a  DE  weapon  of  three-times  the  power  and 
dwell  time  of  the  example  may  very  well  cause  enough  structural  damage  to  “bring  about 
catastrophic  aerodynamic  instability  which  would  tear  the  missile  apart”  [6]. 

The  first  laser,  constructed  by  Theodore  H.  Maiman  in  1960,  was  a  pulsed  ruby 
laser  with  approximately  one  watt  of  output  power  [4].  Almost  immediately  following 
this  initial  invention,  the  U.S.  Department  of  Defense  (DOD)  heavily  invested  into  the 
research  and  development  of  lasers  as  operational  military  weapons  systems.  Several 
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military  entities  have  since  developed  multiple  HEL  weapons  systems  for  operational 
deployment.  The  U.S.  Air  Force  has  constructed  the  Advanced  Tactical  Laser  and  the 
Airborne  Laser  (ABL)  test-bed,  both  of  which  employ  aircraft-mounted  systems  utilizing 
chemical  oxygen  iodine  lasers  (COILs)  with  output  powers  of  100  kW  and  one  megawatt 
(MW),  respectively  [4], 

The  U.S.  Navy  has  also  developed  several  sophisticated  weapons  platforms, 
including  the  Mid-Infrared  Advanced  Chemical  Laser  (MIRACL),  the  Laser  Weapon 
System  (LAWS),  and  the  Maritime  Laser  Demonstrator  (MLD).  The  MIRACL  system, 
first  made  operational  in  1980,  employs  a  deuterium-fluoride  laser  with  ~one  MW  output 
power,  while  LAWS  and  MLD  employ  a  solid-state  fiber  laser  and  a  solid-state  slab  laser 
with  output  power  of  32  kW  and  ~100  kW,  respectively  [4].  The  latter  laser  system,  the 
MLD,  also  became  the  first  to  be  tested  on  and  powered  by  a  ship  in  2011  onboard  the 
U.S.  Navy  defense  test  ship  USS  Paul  Foster  (DD-964).  The  most  significant 
achievement  in  the  implementation  of  a  DE  weapon  system  onboard  a  naval  platform, 
however,  is  the  slated  operational  deployment  of  the  aforementioned  LAWS  on  the  USS 
Ponce  (LPD-15)  in  the  Persian  Gulf  during  summer  2014  [7]. 

B,  ADVANTAGES  OF  HELS 

Conceptually,  directed  energy  weapons  systems  have  significant  superiority  to 
most  conventional  weapons  systems.  The  main  advantages  of  DE  weapons  lie  in  the 
precision  and  speed  of  payload  delivery,  low  cost  per  engagement,  and  an  essentially 
unlimited  magazine  [4].  A  laser  is  capable  of  focusing  to  a  spot  several  centimeters  in 
diameter  on  a  target  from  a  range  of  up  to  10  km  [6].  Unlike  most  conventional  weapons, 
such  as  explosive  missiles  and  kinetic  rounds,  DE  weapon  precision  also  greatly  reduces 
the  possibility  of  collateral  damage.  Furthermore,  due  to  the  speed  of  light  delivery  of 
damage,  a  directed  energy  weapon  has  the  potential  of  engaging  and  destroying  a  target 
soon  after  it  is  detected,  whereas  munitions  from  a  conventional  weapon  must  travel 
toward  the  target  at  a  much  slower  rate  and  thus  engage  it  at  a  close  range.  Radar  may 
also  be  utilized  for  multi-target  tracking,  allowing  for  the  continuous  adjustment  of  the 
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HEL  to  maintain  the  beam  on  a  single  target  while  predieting  the  future  loeation  of  any 
other  targets  for  subsequent  engagement. 

DE  weapons  systems  are  eomprised  of  all-eleetrie  eomponents,  and  eomplete 
integration  into  a  naval  ship’s  power  grid  provides  for  a  magazine  potentially  limited 
only  by  the  ship’s  power.  All-eleetrie  DE  weapons  are  thus  by  design  a  relatively  eost- 
effeetive  platform,  and  several  seeond  “shots”  from  an  HEE  are  only  theoretieally  tied  to 
the  price  per  gallon  of  the  platform’s  fuel.  Although  the  initial  overhead  cost  of  DE 
weapon  design  and  implementation  has  caused  significant  postponement  in  a  budget- 
constrained  environment,  the  cost  per  shot  compared  to  an  advanced  tactical  missile 
makes  the  HEE  a  worthy  long-term  investment. 

C.  HEL  TECHNOLOGY 

The  primary  types  of  DE  weapons  the  U.S.  Navy  is  considering  for  employment 
on  expeditionary  naval  platforms  are  solid-state  lasers  (SSEs)  and  free-electron  lasers 
(EEEs)  [1].  Regardless  of  specific  platform,  the  ideal  HEE  weapon  system  should  exhibit 
several  design  traits.  The  system  should  have  high  wall-plug  efficiency  to  generate  the 
most  optical  power  with  as  little  input  power  as  possible.  The  HEE  should  utilize  a  gain 
medium  with  high  thermal  conductivity  for  rapid  heat  removal  from  the  system.  The 
lasing  medium  should  also  be  resilient  against  heat  and/or  dielectric  breakdown  to  allow 
for  scaling  to  higher  power  without  permanent  damage  to  the  weapon.  Einally,  the  output 
optical  wavelength  should  reside  in  a  region  of  atmospheric  transparency,  so  as  to 
minimize  the  amount  of  attenuation  of  the  laser  energy  and  to  maximize  the  amount  of 
energy  delivered  to  the  target. 

1,  Solid-State  Lasers 

SSL  technology  utilizes  a  solid  substrate  that  holds  the  dopant  in  which  the  lasing 

transitions  occur.  Because  this  technology  utilizes  a  material  medium,  substrate  damage 

threshold  limits  output  power  [4].  SSLs  have  several  advantages  over  EEEs,  including 

their  relatively  compact  and  lightweight  design  (for  low  or  modest  power  outputs), 

relative  power  efficiency,  and  relative  operational  safety  in  the  absence  of  potential  stray 

radiation.  Most  importantly,  SSLs  are  more  technically  mature  than  EEEs,  and  can 
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potentially  see  more  near  term  deployment  availability.  However,  SSLs  also  have  several 
disadvantages  to  FELs,  ineluding  the  laek  of  sealability  to  MW-elass,  a  fixed  optieal 
wavelength  primarily  determined  by  the  gain  medium,  and  diffieulty  in  obtaining 
aeeeptable  beam  quality  for  high-energy  applieations  [1],  [4],  Furthermore,  DE 
applieations  require  the  eombination  of  multiple  solid-state  lasers  to  attain  the  neeessary 
power  output  to  infliet  damage  on  target,  whieh  inereases  their  size  and  weight  and 
degrades  beam  quality. 

The  two  types  of  SSEs  whieh  have  been  employed  as  taetieal  weapons  by  the 
DOD  are  slab  lasers  and  fiber  lasers  [1].  A  slab  laser  eouples  pump  light  with  laser 
energy  through  one  side  of  a  reetangular  medium  slab,  while  removing  heat  from  the 
opposite  side  of  the  slab.  A  fiber  laser,  on  the  other  hand,  utilizes  a  long,  thin,  eylindrieal 
rod  medium,  eovered  by  eladding,  to  eouple  light  to.  While  slab  lasers  have  the 
geometrieal  advantage  of  a  larger  aspeet  ratio  whieh  allows  for  more  light  eoupling  (as 
output  power  is  proportional  to  aspeet  ratio),  fiber  lasers  are  more  lightweight,  power 
effieient,  durable  and  have  eomparable  output  power  to  slab  lasers  with  potentially  better 
beam  quality  in  a  smaller  deviee. 

Furthermore,  optieal  waveguides  in  fiber  lasers  span  many  meters  in  length, 
providing  a  mueh  larger  surfaee  area  that  allows  for  easier  eooling  than  a  slab  laser. 
While  sueh  advantages  indieate  that  the  fiber  SSE  is  a  theoretieally  superior  teehnology 
to  the  slab  SSE,  the  latter  has  aetually  demonstrated  higher  output  power  with  the  ~105 
kW  of  the  MED,  as  opposed  to  the  ~33  kW  fiber  SSE  EAWS  [1].  Ultimately,  however, 
the  most  signifieant  issue  with  SSEs,  whieh  the  FEE  may  be  able  to  address,  lies  in  the 
diffieulty  of  eoherent  beam  shaping  from  eombining  multiple  SSEs,  whieh  is  neeessary  to 
attain  higher  output  power  [4]. 

2,  Free-Electron  Lasers 

Free-eleetron  lasers  utilize  free  (unbound)  eleetrons  to  generate  light  [4].  An 
eleetron  gun,  eonsisting  of  an  eleetron-generating  eathode  and  an  aeeelerator,  ereates  the 
initial  eleetron  beam.  While  several  types  of  eathodes  ean  be  employed  to  generate 
eleetron  beams,  thermionie,  or  “hof  ’  eathodes  are  preferred  for  FEE  applieations  due  to 
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their  much  higher  power  density  [4],  The  electrons  are  accelerated  to  relativistic  speeds 
within  a  linear  accelerator,  and  subsequently  pass  through  an  undulator,  consisting  of  a 
series  of  magnets  which  generate  an  alternating  magnetic  field,  causing  the  electrons  to 
“wiggle”  and  emit  light.  The  interaction  of  the  electrons  with  the  light  in  the  resonator 
produces  stimulated,  coherent  emission.  The  wavelength  of  the  light  produced  depends 
on  the  energy  of  the  electrons  and  the  properties  of  the  undulator.  Furthermore,  wall-plug 
efficiency  of  the  system  can  be  increased  by  recovering  most  of  the  energy  of  spent 
electrons  [4]. 

Because  FELs  utilize  free  electrons  as  the  gain  medium,  elimination  of  waste  heat 
is  not  an  issue  and  thus  their  output  power  is  thought  to  be  scalable  to  megawatt  class. 
Furthermore,  the  FEE  optical  wavelength  is  tunable,  and  thus  can  be  chosen  for  optimal 
atmospheric  transparency.  An  EEE  can  also  obtain  excellent  beam  quality  at  high  output 
power  for  long-distance  propagation  and  a  relatively  energy-efficient  design.  However, 
EEEs  tend  to  be  rather  large  in  size  ( ~  20  meters  long),  heavy,  and  expensive  to  design 
[1].  EEEs  also  require  shielding  to  contain  radiation  produced  by  stray  electrons,  and  the 
overall  technology  is  relatively  less  mature  compared  to  SSEs. 

D.  HEL  MISSIONS  AND  PLATFORMS 

The  missions  of  a  DE  weapon  are  largely  dependent  on  its  output  power  level. 
Directed  energy  weapons  with  lower  power  output  (10-100  kW),  such  as  the  SSE,  are 
better  suited  against  slower  moving  targets,  such  as  small  watercraft  and  aircraft.  Such 
DE  weapons  could  immobilize  boats  and  down  aerial  drones,  as  well  as  offer  the 
potential  of  “soft  kills”  by  blinding  incoming  missiles.  Megawatt-class  DE  weapons,  such 
as  the  EEE,  will  be  primarily  slated  for  cruise  missile  defense.  The  current  EEE  DOD 
project,  the  Innovative  Naval  Prototype,  maintains  a  planned  100  kW  of  output  power 
and  is  a  stepping  stone  in  the  future  transition  to  the  MW-class  devices. 

The  HEE  weapons  systems  of  both  kW  and  MW-class  could  be  fittingly 
employed  on  several  naval  platforms,  such  as  the  future  all-electric  ships,  littoral  combat 
ships  (ECS),  guided-missile  destroyers  (DDG-51  and  DDG-1000),  and  the  future  Eord 
class  aircraft  carriers,  such  as  the  Gerard  R.  Ford  (CVN-78).  Due  to  weight  and  size 
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constraints,  only  SSLs  could  be  potentially  employed  on  Naval  aviation  platforms  sueh  as 
F/A-18E/F  Hornet,  F-35  Joint  Strike  Fighter  (JSF)  [8],  and  Marine  Corps  helieopter 
platforms  sueh  as  the  UH-IY  Huey  and  AH-IZ  Cobra  [6]. 
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III.  ATMOSPHERIC  PROPAGATION  OVERVIEW 


The  effectiveness  of  HEL  weapons  systems  can  be  greatly  affected  by  the 
medium  through  which  the  laser  energy  propagates.  Because  this  medium  is  the  Earth’s 
atmosphere,  one  needs  to  understand  the  various  optical  properties  and  physical 
principles  associated  with  this  environment  to  be  able  to  develop  systems  that  can  operate 
effectively  within  it. 

A,  ATMOSPHERIC  EXTINCTION 

Molecules  and  aerosols  suspended  in  the  atmosphere  absorb  and  scatter 
propagating  laser  energy  as  it  travels  to  reach  the  target  [4].  Absorption  occurs  when 
photons  from  incident  radiation  (such  as  laser  light)  transfer  their  energy  to  atmospheric 
constituents,  thereby  heating  the  local  atmosphere.  Scattering  occurs  when  incident 
photons  are  redirected  from  the  beam  in  random  directions  by  atmospheric  particulates. 

These  combined  effects  are  illustrated  by  Beer’s  Eaw  in  Equation  (1),  which 
describes  the  transmission  of  light  through  a  linear  medium,  following 

/■(z)  =  (1) 

In  Equation  (1),  Prepresents  the  total  power  delivered  to  a  plane  located  a 
distance  z  from  the  source  and  is  the  initial  output  power  at  the  beam  director.  The 

variable  e  represents  the  total  extinction  coefficient  due  to  atmospheric  absorption  and 
scattering  (in  units  of  inverse  length),  and  is  further  described  by  Equation  (2)  as 

e  =  CC^+CCa^  Pa-  (2) 

In  Equation  (2),  a  refers  to  the  absorption  coefficients  and  refers  to  the 
scattering  coefficients,  while  subscripts  m  and  a  denote  the  molecular  and  aerosol 
contributions,  respectively.  An  increase  in  e  results  in  greater  attenuation  rate  of  laser 
light  by  the  atmosphere. 

1.  Molecular  Effects 

Table  1  displays  a  list  of  the  major  molecular  constituents  which  make  up  the 
Earth’s  atmosphere. 
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The  Composition  of  the  Atmosphere" 


Permanent  constituents  Variable  constituents 


Constituent 

%  by  volume 

Constituent 

%  by  volume 

Nitrogen  (N2) 

78.084 

Water  vapor  (H2O) 

0-0.04 

Oxygen  (O2) 

20.948 

Ozone  (O3) 

0-12  X  10-* 

Argon  (Ar) 

0.934 

Sulfur  dioxide  (SO2)* 

0.001  X  10“-^ 

Carbon  dioxide  (CO2) 

0.036 

Nitrogen  dioxide  (NO2)* 

0.001  X  10--* 

Neon  (Ne) 

18.18  X  lO--* 

Ammonia  (NH3)'' 

0.004  X  10-* 

Helium  (He) 

5.24  X  10-* 

Nitric  oxide  (NO)* 

0.0005  X  10-* 

Krypton  (Kr) 

1.14  X  10-^ 

Hydrogen  sulfide  (H2S)* 

0.00005  X  10- 

Xenon  (Xe) 

0.089  X  10-^ 

Nitric  acid  vapor  (HNO3) 

Trace 

Hydrogen  (H2) 

0.5  X  lO-* 

Chlorofluorocarbons 

Trace 

Methane  (CH4) 

1.7  X  10-^ 

(CFCI3,  CF2CI2 

Nitrous  oxide  (N2O)* 

0.3  X  lO--’ 

CH3CCI3,  ecu,  etc.) 

Carbon  monoxide  (CO)* 

0.08  X  10-“* 

“After  the  U.S.  Standard  Atmosphere  (1976)  with  modifications. 
*  Concentration  near  the  earth’s  surface. 


Table  1.  Molecular  composition  of  Earth’s  atmosphere.  From  [9]. 


Ambient  air  is  predominately  composed  of  N2  and  O2,  both  of  which  are  linear, 
homonuclear  molecules  with  symmetric  molecular  bonds  between  each  of  their  two 
atoms.  Because  these  molecules  only  possess  one  normal  vibrational  mode  that  does  not 
alter  their  overall  dipole  moments,  propagating  optical  energy  does  not  cause  these 
molecules  to  vibrate.  Thus,  N2  and  O2  also  do  not  absorb  electromagnetic  radiation  in  the 
infrared  spectral  region  [9]. 

Instead,  other  atmospheric  constituents,  such  as  H2O  vapor  and  CO2,  dominate  the 
infrared  atmospheric  absorption  spectrum.  Even  though  these  molecules  contribute  to  a 
very  small  percentage  of  the  total  atmospheric  composition,  they  are  asymmetrically- 
shaped  and  exhibit  multiple  vibrational  modes,  making  them  strong  absorbers  of 
electromagnetic  radiation  [9].  The  molecular  absorption  is  extremely  intricate  and  ideal 
laser  performance  requires  that  the  spectral  width  of  the  beam  fits  within  a  wavelength 
window  containing  minimum  absorption  features.  A  typical  molecular  absorption 
spectrum  for  a  wavelength  range  of  0.5  pm  to  5.0  pm  is  shown  in  Figure  1.  The  gaps  in 
the  data  at  ranges  ~2.5  to  3  pm,  ~4.25  to  4.5  pm  and  -4.75  to  5.0  pm  are  due  to  the 
absorption  coefficient  nearing  infinity,  which  represents  wavelength  ranges  of  high 
absorption. 
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Molecular  Absorption  Spectrum 


Figure  1.  Typical  molecular  absorption  spectrum  for  the  atmosphere.  Gaps  in 
the  data  correspond  to  high  molecular  absorption  coefficient  . 


Figure  2  displays  a  typical  molecular  scattering  spectrum  for  the  atmosphere  over 
a  wavelength  range  of  0.5  pm  to  5.0  pm.  Molecular  scattering  is  described  by  an  optical 
phenomenon  known  as  Rayleigh  scattering,  in  which  shorter  wavelengths  of 
electromagnetic  energy  are  scattered  more  readily  than  longer  wavelengths.  Rayleigh 
scattering  will  be  further  discussed  in  Chapter  IV,  Section  B.4. 


Figure  2.  Typical  molecular  scattering  spectrum  for  the  atmosphere.  As  a 
result  of  Rayleigh  scattering  (discussed  in  Chapter  IV),  shorter 
wavelengths  of  light  scatter  more  readily  than  longer  wavelengths. 
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Most  of  the  photons  scattered  by  molecules  have  the  same  wavelength  as  the 
incident  photons;  only  their  direction  of  propagation  changes.  This  principle  readily 
describes  the  phenomenon  of  the  blue  daytime  sky  and  red  sunsets.  Because  blue  light 
more  readily  undergoes  molecular  scattering  than  red  light,  the  daytime  sky  appears  blue 
[9].  At  sunset,  the  path  length  of  light  from  the  setting  sun  to  the  viewer  through  the 
atmosphere  is  longer,  and  by  the  time  the  light  reaches  the  viewer’s  eyes  most  of 
the  shorter  wavelength  light  has  been  scattered,  leaving  only  the  longer  wavelength  red 
light  [9]. 

2,  Aerosol  Effects 

Absorption  and  scattering  effects  on  light  by  aerosols  are  similar  in  nature  to 
those  by  molecules.  However,  the  extent  of  these  effects  depends  greatly  on  aerosol  size, 
shape  and  refractive  index,  the  latter  primarily  defined  by  the  material  properties  of  the 
individual  particles  (water,  dust,  etc.).  Figure  3  displays  a  plot  of  the  aerosol  absorption 
spectrum,  while  Figure  4  displays  a  plot  of  the  aerosol  scattering  spectrum,  both  for  a 
0.5  pm  to  5.0  pm  wavelength  range. 

Aerosol  Absorption  Spectrum 
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Figure  3.  Typical  aerosol  absorption  spectrum  for  the  atmosphere.  The  peak  in 
aerosol  absorption  at  ~  3  pm  corresponds  to  presence  of  particulates 
which  specifically  absorb  at  those  wavelengths. 
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Aerosol  Scattering  Spectrum 
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Figure  4. 


Typical  aerosol  scattering  spectrum  for  the  atmosphere.  Aerosols 
scatter  most  wavelengths  relatively  equally. 


Because  most  aerosols  are  similar  in  size  to  the  wavelength  of  infrared  or  visible 
light,  they  absorb  most  wavelengths  approximately  equally  [5].  The  notable  increase  in 
scattering  at  ~  3  pm  may  be  due  to  particulates  such  as  sulfate  (SO2)  or  sea  salt  aerosols, 
which  may  absorb  more  at  this  specific  wavelength  because  of  their  size. 

Similarly,  aerosol  particulates  also  scatter  most  wavelengths  in  the  infrared  and 
visible  range  relatively  equally,  which  explains  why  clouds  and  haze  appear 
approximately  white  to  the  human  eye  [9]. 

B,  ATMOSPHERIC  TURBULENCE 

Turbulence  is  a  term  generally  describing  a  fluid  flow  regime  within  a  medium  in 
which  the  flow  exhibits  an  irregular,  chaotic  and  dissipative  pattern,  mostly  in  the  form  of 
randomly-varying  eddies  [10].  In  the  atmosphere,  the  turbulent  flow  of  air  is  considered 
to  be  small-scale  and  primarily  the  result  of  convective  air  motion,  vertical  temperature 
differences,  and  wind  shear. 

The  major  effects  of  turbulence  on  electromagnetic  energy  propagation  stem  from 
the  spatial  and  temporal  fluctuations  in  the  refractive  index  of  air  [10].  The  relationship 
between  absolute  temperature  T  (in  K),  wavelength  X  (in  pm),  pressure  P  (in  mbar) 

and  the  refractive  index  of  air  n  can  be  described  by  [10] 
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Sufficient  turbulence  in  the  local  atmosphere  leads  to  variations  in  the  local 
ambient  temperature,  which  in  turn  leads  to  variations  in  the  refractive  index.  For  HEL 
weapons  systems  to  function  most  effectively,  the  propagation  of  laser  energy  through 
the  atmosphere  requires  that  the  light  rays  emitted  by  the  laser  maintain  transverse 
coherence  so  that  they  constructively  interfere  to  form  a  high  irradiance  spot  size  on 
target.  However,  refractive  index  variations  due  to  turbulence  cause  different  parts  of  the 
laser  beam  to  experience  different  path  lengths  as  the  beam  propagates,  which  disrupts 
constructive  interference  at  the  target  and  can  even  cause  the  beam  to  spread  out  and/or 
scintillate  (i.e.,  breakup  into  beamlets)  [10]. 


The  Fried  parameter  (units  of  length)  can  be  used  to  estimate  the  effects  of 
turbulence  on  laser  propagation.  The  parameter  defines  a  circular  diameter  over  which 
the  laser  beam  maintains  coherence  (no  significant  change  in  phase)  in  the  transverse 
direction  throughout  the  propagation  distance  [4].  In  other  words,  estimates  the 
maximum  transverse  separation  at  the  source  between  two  individual  light  rays  within 
which  the  two  rays  will  constructively  interfere  at  the  target.  A  lower  value  implies 

stronger  turbulence.  If  the  value  is  significantly  smaller  than  the  beam  director  size, 
the  laser  beam  will  break  up  into  many  smaller,  incoherently  radiating  beamlets. 


The  relationship  between  ,  the  refractive  structure  constant  C] ,  the  wavelength 
of  light  /I,  and  the  distance  to  target  t  can  be  quantified  by 


Tq  »  0.33 


/5(y2 


(4) 


For  typical  wavelength  and  target  ranges  for  HEF  applications,  can  range  from  ~one 
centimeter  in  strong  turbulence  to  lO’s  of  centimeters  in  weak  turbulence.  The  Cl  value 
(in  units  of  m'^^^)  characterizes  the  strength  of  the  turbulence.  The  Cl  value  was  first 
theoretically  parameterized  by  Audrey  Kolmogorov  and  has  been  measured  at  numerous 
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locations  worldwide  [4],  The  signifleance  of  the  Cl  value  stems  from  Kolmogorov’s 
original  theory  for  turbulent  flow,  whieh  has  formed  the  basis  for  all  eontemporary 
turbulenee  theories  [4], 

The  basie  idea  behind  Kolmogorov’s  theory  is  this:  Turbulent  fluid  flow  is 
governed  by  the  Navier-Stokes  equations,  but  signifieant  ehallenges  he  in  solving  these 
equations  within  a  system  as  expansive  as  the  atmosphere.  Kolmogorov’s  theory 
therefore  utilizes  a  statistieal  approach  in  describing  the  flow  of  kinetie  energy  from 
large-scale  eddies,  whieh  can  be  tens  to  hundreds  of  meters  in  size,  into  smaller  eddies, 
whieh  are  typieally  millimeters  to  eentimeters  in  size  [10].  Eddies  within  this  range  of 
sizes,  whieh  is  defined  as  the  inertial  subrange,  are  assumed  to  be  relatively 
homogeneous  and  isotropie  within  smaller  regions  of  spaee  [4],  and  thus  properties  sueh 
as  the  refraetive  index  vary  over  spaee  and  time  in  identieally-distributed  inerements. 
This  assumption  allowed  for  Kolmogorov  to  utilize  a  strueture  funetion  in  deseribing 
turbulence,  whieh  is  represented  by  the  Cl  value. 

In  atmospheric  modeling,  an  increase  in  the  Cl  value  yields  an  increase  in  the 

amount  of  atmospheric  turbulence.  Furthermore,  the  parameter  tends  to  be  a  strong 
function  of  altitude,  and  generally  deereases  at  higher  altitudes.  Weak  atmospherie 
turbulenee  would  exhibit  a  value  of  C]  «10"'’m“^'%  while  strong  turbulence  would 
oceurif  Cl  «10-'^m-"'\ 

C.  THERMAL  BLOOMING 

Thermal  blooming  deseribes  a  non-linear  phenomenon  in  whieh  propagating  laser 
light  deposits  a  small  portion  of  its  energy  along  the  beam  path  into  the  local  atmosphere. 
As  the  laser  energy  is  absorbed  by  the  molecules  and  aerosols  in  the  air,  the  ambient  air 
temperature  increases,  creating  localized  gradients  in  air  density,  subsequently  altering 
the  refraetive  index  of  the  medium  [11]  (see  Equation  3).  A  suffieient  ehange  in  the 
refraetive  index  effeetively  ereates  a  diverging,  eoneave  lens  within  the  volume  of  the 
laser  beam,  whieh  exhibits  a  greater  lens  thiekness  and  thus  optieal  density  at  the  edges 
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than  at  the  center  [5].  Ultimately,  this  lens  causes  the  beam  to  “bloom,”  or  diverge  from 
its  propagation  path,  effectively  reducing  the  irradiance  on  target. 

1.  Paraxial  Wave  Equation 

The  effects  of  thermal  blooming  on  HEL  beam  propagation  can  be  estimated  by 
first  assuming  that  the  electric  field  of  the  propagating  electromagnetic  energy  takes  the 
form  of 

E{r)  =  y/{r)e^’'\  (5) 

where  is  the  slowly-varying  complex  envelope  and  the  harmonic  factor  has 

been  suppressed.  Here,  k-lnik  is  the  wave  number,  X  is  the  laser  wavelength,  z  is 
the  primary  direction  of  propagation,  and  the  real  part  of  E  specifies  the  electric  field. 
The  irradiance  of  the  energy  in  SI  units  thus  takes  the  form  of 

/(r)  =  ^£-£  =  ^rV,  (6) 

where  y/*  represents  the  complex  conjugate  of  wave  function  y/  .  Equation  (5)  represents 
a  solution  to  the  paraxial  wave  equation,  which  follows  as  [11] 

^  =  -^y/  +  ^V\y/  +  {ik5n)y/ ,  (7) 

where  -ey/ !  2  represents  atmospheric  extinction  (see  section  III  A.),  iy\y/ /  2k 
represents  transverse  diffraction,  V\=  !  dx^ !  dy^  is  the  transverse  Eaplace 

operator,  and  the  {ik5n)y/  component  describes  the  effects  of  varying  the  refractive 
index  5n  due  to  thermal  blooming  (and  in  principle,  turbulence). 

As  stated  previously,  a  change  in  air  temperature  has  a  direct  effect  on  the 
refractive  index  of  the  absorbing  medium,  described  by 

5n  =  —^T  =  -'^^^T ,  (8) 

dT 

where  is  the  averaged  refractive  index  of  air,  is  the  ambient  temperature,  and  AT 
represents  the  increase  in  temperature  due  to  absorption. 
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For  most  HEL  weapon  systems,  it  is  common  to  assume  that  the  rate  of  heating  is 
too  small  to  create  pressure  differences  (assumed  to  be  nearly  instantaneous)  [11], 
resulting  in  constant-pressure,  or  isobaric,  heating.  The  temperature  rise  due  to  laser 
heating  under  these  conditions  can  be  described  using  the  energy  balance  equation,  given 
by  [11] 


^^T 


a 


■/-(v.v^)Ar+ 


K 


pC 


■V^ 


1 

p  J 


AT. 


(9) 


dt  pCp 

In  Equation  (9),  the  a  I  /  pC^  term  represents  laser  heating,  where  a  is  the 
atmospheric  absorption  coefficient,  p  is  the  air  density,  is  the  constant-pressure 
specific  heat,  and  /  is  the  laser  irradiance.  The  (v»V^)Arterm  represents  the  effects  of 
energy  exchange  by  convection,  where  v  is  the  wind  velocity,  is  the  transverse 
gradient  operator.  The  /  pC^lSJ  term  represents  conduction  of  heat  through  the 

air,  where  K  is  the  thermal  conductivity. 


In  most  HEE  propagation  scenarios,  convection  dominates  conduction  except  for 
conditions  involving  absolutely  calm  winds.  If  one  assumes  that  the  laser  beam’s  pulse 
repetition  frequency  (PRE)  remains  low  enough  to  allow  the  locally  heated  air  to 
be  removed  from  the  beam  path  by  either  wind  or  beam  motion,  then  the  conduction 
term  can  essentially  be  ignored  in  what  is  defined  as  convection-dominated  thermal 
blooming  [11]. 


2.  Steady-State  Phase  Distortion  without  Wind 


In  the  case  of  a  stationary  beam  and  no  wind  (|v  =  O) ,  the  thermal  blooming  due 


to  absorption  approaches  the  steady-state  conduction-dominated  scenario.  The  magnitude 
of  conduction-dominated  thermal  blooming  on  the  laser  beam  can  be  evaluated  using  the 
dimensionless  scaling  parameter  given  by  [12] 


rij^Paz  . 

D..= - - - ^ ,  az  <Si\. 

^  InKna 


(10) 
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where  from  Equation  (8)  we  obtain  that  nj.=dn/dT  =  -(^n^-T^)/ T^,  P  is  the  laser 

power,  z  is  the  distanee  that  the  beam  has  propagated  through  the  absorbing  medium, 
and  a  is  the  beam  radius.  The  on-axis  decrease  in  beam  intensity  due  to  conduction- 
dominated  thermal  blooming  can  be  approximated  by  [11]. 

For  the  contribution  of  conduction  to  thermal  blooming  to  be  ignored,  the  wind 
speed  |v|  must  satisfy  the  condition  where  x  =  KI  pC^  represents  the 

thermal  diffusivity  of  the  medium  (z  ~  0.2cm^  /  s  for  air  at  300  K)  [11].  Because  the 
ratio  of  the  thermal  diffusivity  to  the  beam  radius  is  virtually  negligible  in  comparison  to 
wind  speed  for  all  scenarios  but  those  involving  absolutely  no  wind,  this  condition  is 
nearly  always  satisfied  for  laser  propagation. 

3.  Steady-State  Phase  Distortion  with  Wind 

Phase  distortion  due  to  thermal  blooming  in  the  presence  of  wind  has  a  very 
specific  effect  on  the  laser  beam  irradiance  on  target,  illustrated  by  Figure  5. 


incident  refractive  distorted 

INTENSITY  TEMPERATURE  INDEX  INTENSITY 


Figure  5.  The  effects  of  steady-state  thermal  blooming  with  a  transverse  wind, 

from  [11]. 
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As  the  transverse  wind  eomponent  moves  aeross  the  beam  propagation  path,  it 
eauses  the  loeal  air  temperature  to  inerease  in  the  downwind  direetion.  As  temperature 
change  is  inversely  proportional  to  both  the  change  in  air  density  and  the  refractive  index, 
both  quantities  decrease  in  the  same  direction  as  the  wind.  Because  light  bends  in  the 
direction  of  an  increasing  refractive  index,  the  laser  beam  shifts  its  peak  intensity  in  the 
upwind  direction. 

This  effect  can  be  roughly  quantified  by  assuming  steady-state  thermal  blooming, 
in  which  5Ar  /  =  0  in  Equation  (9),  and  that  the  wind  velocity  is  entirely  along  the 


transverse  x-direction,  so  \-vi.  The  phase  shift  (in  rad)  from  thermal  blooming  can  then 
be  written  as  [11] 

(Pg[x,y,z)  =  ik{Sndz  =  -— — -  \ldx' ,  (11) 

0  PCpVi 

where  the  beam  is  assumed  to  be  collimated  along  the  z-direction.  If  we  further  assume 
an  irradiance  distribution  for  a  Gaussian  beam.  Equation  (11)  can  be  solved  analytically 
to  obtain  [11] 


g)p^{x,y,z) 


exp  -[^y  I  of  [l-l-e//(x/a)] , 
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where  the  maximum  phase  shift  across  the  beam  Acp^  can  be  expressed  as 

(  1  ^ 

~  ^  j— 
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The  parameter  represents  the  Bradley-Hermann  [13]  distortion  number  and  is 
utilized  as  a  quantitative  measure  for  the  intensity  of  thermal  blooming  [11],  Blooming 
beeomes  a  significant  complication  to  the  propagation  of  laser  energy  when  >  25  [4], 

The  effects  of  phase  distortion  on  the  laser  beam  can  be  specifically  inferred  from 
the  series  expansion  in  Equation  (13).  The  leading  (x/  a)  term  produces  the  deflection  of 

the  laser  beam  along  the  wind  direction  in  the  x-direction  and  the  leading  (y  /  a)^  term 
produces  the  beam  divergence  in  the  y-direction  [11], 


4.  Laser  Beam  Irradiance  Pattern  Due  to  Thermal  Blooming 


The  thermal  blooming  irradiance  pattern  can  be  estimated  for  the  phase 
accumulation  by  solving  the  paraxial  wave  equation  in  the  ray-optics  limit  using 
perturbation  theory.  To  the  first-order  perturbation,  the  solution  follows  as  [11] 


I{x,y,z) 

I 


=  exp 


(16) 


where  is  the  unperturbed,  Gaussian,  collimated  beam  profile  and  is  the  collimated 
beam  irradiance  distortion  parameter  with  the  inclusion  of  wind  effects,  taking  the  form 

of[ll] 


^  _{n,-\)aPz^ 
nn^pC^Tva^ 


(17) 


In  general,  an  increase  in  the  beam  distortion  parameter  represents  a  greater  effect 
of  thermal  blooming  on  the  laser  beam  and  results  in  an  irradiance  pattern  taking  the  form 
of  a  crescent  shape,  which  is  illustrated  in  Figure  6. 
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Figure  6.  The  general  transverse  laser  beam  profile  due  to  the  effects  of 
steady-state  thermal  blooming.  The  presence  of  wind  moving  from  left 
to  right  in  the  figure  results  in  the  typical  crescent  shape  shown.  The 
color  map  indicates  the  irradiance  of  the  laser  beam,  with  the  dark  red 
color  representing  the  highest  irradiance  and  the  dark  blue  the  lowest. 
The  scale  of  beam  deflection  depends  upon  the  distance  to  the  target, 
beam  power,  atmospheric  extinction,  and  other  parameters  as  described 

in  the  text. 
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IV.  ATMOSPHERIC  CODE  OVERVIEW 


Effective  employment  of  HEL  weapon  systems  requires  knowledge  of  the 
atmospheric  optical  properties  and  how  these  properties  affect  laser  propagation.  The 
atmospheric  parameters  that  determine  these  optical  properties  are  extremely  complicated 
and  depend  on  pressure,  temperature,  humidity,  aerosol  size  composition,  etc. 
Eurthermore,  all  of  these  parameters  are  functions  of  elevation,  time  of  day,  and  season. 

Due  to  these  complications,  it  is  difficult  to  fully  map  out  this  multi-dimensional 
parameter  space  for  varying  altitudes  with  experimental  data.  Thus,  a  means  of 
generating  optical  properties  based  upon  known  physical  principles  for  a  large  range  of 
conditions  must  be  employed,  which  can  then  be  used  to  evaluate  laser  beam 
propagation. 

A.  MOLECULAR  ABSORPTION: 

Among  the  primary  atmospheric  optical  properties  that  govern  laser  beam 
propagation  are  the  four  components  of  the  extinction  coefficient  (See  Chapter  II,  Section 
A).  A  combination  of  measurements  and  quantum  theory  is  utilized  to  determine  the  first 
component,  the  molecular  absorption  coefficient  ,  which  represents  the  absorption  of 
laser  energy  by  molecules  in  the  atmosphere. 

Within  the  electromagnetic  spectrum,  the  position  and  strength  of  absorption 
features  at  each  wavelength  arise  from  quantized  molecular  energy  levels.  The 
atmospheric  codes  which  are  employed  in  atmospheric  modeling  query  a  database  of 
absorption  features  (such  as  the  HITRAN  database),  modify  the  line  shapes  based  upon 
temperature  and  pressure,  and  then  add  them  together  to  get  the  absorption  coefficient  at 
a  specific  wavelength.  This  method  is  based  on  assumptions  regarding  the  natural  line 
shape  and  how  temperature  and  pressure  affect  the  shape. 
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1.  Assumptions 

The  general  form  of  an  absorption  ooeffleient  line  shape  is  represented  by 

a,=NSf(v-v^),  (18) 

where  N  represents  the  number  of  atoms  per  unit  volume,  S  is  the  line  strength, 
/  (v  -  Vg )  is  the  normalized  line  shape  at  frequeney  v  ,  and  is  the  central  frequency  of 

the  line,  where  /(v-v^)  dv  —  One  of  the  most  fundamental  assumptions  regarding 

molecular  absorption  modeling  is  that  the  natural  line  shape  for  each  absorption  feature  is 
based  upon  modeling  a  molecular  electron  as  a  damped,  driven  harmonic  oscillator,  with 
the  equation  of  motion  given  by  [9]; 

F  eE  ■■ 

—  =  —  =  r  +  4;7r  r+<^r.  (19) 

In  Equation  (19),  F  represents  the  force  vector,  is  the  mass  of  the  electron,  e 

is  the  charge  of  the  electron,  E  is  the  total  electric  field  (including  the  contribution  from 
the  light),  f  is  the  electron  acceleration,  T ^  is  the  effective  natural  damping  constant,  f 
is  the  electron  velocity,  ^  is  related  to  the  strength  of  the  restoring  force  (which  is  related 
to  the  effective  spring  constant),  and  f  is  the  electron  position.  While  the  potential 
energy  of  molecular  bonds,  also  known  as  the  Morse  potential,  is  not  the  same  as  a 
harmonic  oscillator,  at  lower  energy  levels  it  can  be  reasonably  approximated  by  the 
harmonic  potential,  as  shown  in  Figure  7,  thus  providing  motivation  for  using  Equation 
(19).  This  approximation  is  nominally  satisfied  at  typical  atmospheric  conditions,  as  most 
of  molecular  electrons  remain  in  their  ground  energy  state  at  standard  temperature  and 
pressure  (STP). 
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Figure  7.  Harmonic  oscillator  approximation.  At  low  energy  levels  and  small 
intemuclear  separation,  the  Morse  potential  can  be  adequately 
approximated  by  the  harmonic  potential,  after  [15]. 


Equation  (19)  can  be  used  to  determine  the  complex  refractive  index  [9].  The 
imaginary  part  of  the  refractive  index  is  responsible  for  absorption  and  determines  the 
functional  shape  of  /  ( v  -  )  as 


f{y-K)  = 


1 


(20) 


where  is  the  same  effective  damping  constant  in  Equation  (19).  Equation  (20)  is  in 
the  form  of  a  Lorentzian,  with  E^  also  specifying  the  full-width  half-maximum  of  the 
response,  as  visualized  in  Figure  8. 
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Figure  8.  The  dependence  of  spectral  line  width  /(v  -  and  height  on 
resonant  frequency  and  damping,  respectively,  after  [16]. 


2.  Absorption  Line  Broadening 

The  absorption  spectrum  of  a  single  line  is  broken  down  into  the  shape  and 
strength  of  the  line  as  described  by  Equation  (18).  However,  various  mechanisms  affect 
these  spectral  line  characteristics,  which  are  discussed  below. 


a.  Natural  Line  Width 


The  natural  line  width  of  an  absorption  feature  is  characterized  by  the  damping  of 
the  harmonic  oscillator  vibrations  due  to  energy  emission  of  the  molecule  [9].  The 
Heisenberg  uncertainty  principle  can  be  used  to  estimate  the  width  of  the  absorption  line 
and  thus  provides  a  means  of  estimating  .  It  states 


h 
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(22) 


where  AE  is  the  change  in  energy  of  the  molecule  for  a  particular  transition.  At  is  the 
lifetime  of  the  excited  state,  h  is  the  Planck  constant,  and  Av  is  the  width  of  the  line 
shape. 
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In  the  thermal  infrared  regime,  the  lifetimes  of  isolated  moleeules  are 
At  «  0.1  - 10  s,  in  which  case  the  width  of  the  lineshape  A  v  ~  0.3  Hz,  which  means  the 
damping  constant  T ^  is  similarly  small.  Because  the  natural  line  width  determined  from 

the  uncertainty  principle  and  modeled  by  quantum  mechanics  is  very  small,  the  effects  of 
the  two  broadening  mechanisms  discussed  below  dominate  the  spectral  line  shapes  and 
can  be  neglected. 


b.  Pressure  Broadening 


As  stated  earlier,  two  broadening  mechanisms  dominate  the  characteristic  shape 
of  spectral  lines.  Pressure  broadening  is  a  mechanism  by  which  molecular  collisions  in 
the  atmosphere  transfer  energy  between  molecules,  thereby  reducing  the  lifetime  of 
higher  energy  states  and  increasing  the  effective  damping  of  the  oscillator.  It  is  the 
dominant  broadening  mechanism  in  the  lower  atmosphere  below  ~  20  km  [14]. 
Increased  pressure  causes  more  collisions  that  broaden  the  Lorentzian  absorption  lines, 
taking  the  form  of 


(- 


Tp!  TT 
-Vo)'+r, 
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(23) 


where  is  the  effective  damping  constant  due  to  collisions  and  is  many  orders  of 


magnitude  larger  than  F^  [14].  The  damping  constant  F^  is  determined  from  the  mean 


time  between  inelastic  molecular  collisions  ,  so  that  F^  =  1  /  .  It  also  follows  the 

trend 
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where  F^  represents  the  line  width  at  a  reference  temperature  and  pressure  ,  and 
n  « 1/2  is  determined  empirically  [9]. 


Figure  9,  below  illustrates  the  Lorentz  line  shape  for  three  pressure  values.  For 
typical  vibration-rotation  bands,  a  line  width  of  1.5x10^  Hz  at  a  pressure  of  one  bar  is 
typical  [17]. 
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Figure  9.  Lorentz  line  shape  eaused  by  pressure  broadening,  with  variations  in 
damping  constant  and  atmospheric  pressure,  after  [14], 


Doppler  Broadening 


c. 


Doppler  shifts  also  serve  to  broaden  spectral  line  shapes.  As  laser  energy 
propagates  through  the  atmosphere,  the  light  encounters  molecules  with  a  random 
distribution  of  velocities.  A  Doppler  frequency  shift  due  to  molecular  motion  results  in  a 
further  broadening  of  the  emission/absorption  spectrum,  but  in  a  Gaussian  profile 
described  by  [18] 


(25) 


(26) 


where  is  the  Doppler  line  half-width,  kg  is  the  Boltzmann  constant,  m  is  the  mass  of 
the  molecule  and  c  is  the  speed  of  light. 

Since  a  change  in  temperature  implies  a  change  in  the  velocity  of  molecular 
motion,  an  increase  in  temperature  generally  correlates  to  an  increase  in  Doppler 
broadening. 


28 


d.  Voigt  Profile 

Atmospheric  codes  use  Voigt  profiles  to  combine  the  effects  of  pressure 
(Lorentzian)  and  Doppler  (Gaussian)  broadening  mechanisms.  The  Voigt  line  shape  is 
quantitatively  described  by  [14] 


r  7 


-dy. 


^^'^ToL{x-y)+T 

where  T  =  and  x  =  {y  -Vq)  !  T  The  Voigt  profile  is  illustrated  in  Figure  10. 
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Figure  10.  Collision,  Doppler  and  Voigt  spectral  line  profiles,  after  [14]. 


In  the  shoulders  of  the  line  shape,  the  Voigt  profile  exhibits  a  Lorentzian  shape,  while  a 
resemblance  to  the  Gaussian  shape  is  seen  near  the  center  of  the  line  shape. 

3,  Continuum  Absorption 

Although  the  Voigt  profile  is  a  well-known  representation  of  spectral  line  shapes, 
it  is  not  completely  devoid  of  error.  Continuum  absorption,  a  phenomenon  discovered  in 
1918,  represents  the  discrepancy  between  such  data  obtained  from  a  mathematical 
summation  of  Voigt  lines  and  actual  measurement  data  of  the  molecular  absorption 
spectrum. 
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The  physical  principles  behind  continuum  absorption  remain  debated;  however,  a 
couple  of  theories  attempt  to  reconcile  the  difference  between  experimental  and  modeling 
data.  One  possible  explanation  lies  in  the  hypothesis  that  actual  absorption  line  shapes  are 
not  perfectly  modeled  by  the  Voigt  profile,  and  a  summation  over  thousands  of  such  lines 
gradually  yields  a  rather  significant  margin  of  error.  Another  possible  explanation  lies  in 
the  tendency  of  water  to  form  dimers  in  the  atmosphere.  A  dimer  consists  of  two  H2O 

molecules  loosely-bound  with  a  hydrogen  bond,  exhibiting  measured  dissociation  energy 
of  ~  13.2±0.12kJ / mol  [19].  Incident  laser  light  upon  such  dimers  could  cause  the  beam 
to  deposit  extra  energy  into  the  atmosphere,  resulting  in  an  absorption  spectrum  which  is 
not  properly  modeled  by  the  Voigt  profile. 

Whatever  the  mechanism,  continuum  absorption  varies  slowly  with  wavelength, 
and  thus  atmospheric  codes  often  add  the  contribution  from  continuum  absorption  using 
the  semi-empirical  Clough-Kneizys-Davies  (CKD)  continuum  model  [20]. 

B.  MIE  SCATTERING  THEORY:  «„,/?„,  AND  /?„ 

While  a  combination  of  measurements  and  quantum  theory  is  utilized  to 
determine  the  first  component  of  the  extinction  coefficient,  a  separate  solution  called  Mie 
scattering  is  used  to  determine  the  latter  three  components. 

1.  Theory  Introduction 

Mie  scattering  theory  describes  the  interaction  between  electromagnetic  radiation 
and  small  particles  suspended  in  the  atmosphere,  such  as  dust  and  aerosols  [9].  The 
theory  is  utilized  to  determine  the  molecular  scattering  coefficient  as  well  as  the 

aerosol  absorption  coefficient  and  aerosol  scattering  coefficient  .  Mie  theory 

contains  three  key  assumptions,  the  first  of  which  states  that  all  atmospheric  particles  are 
spherical  in  shape.  The  second  assumption  states  that  the  scattering  particles  are 
homogeneous  in  material  and  thus  can  be  represented  by  a  single  complex  refractive 
index  h.  The  last  assumption  holds  that  the  light  incident  upon  a  particle,  in  this  case 
from  the  laser,  is  in  the  form  of  a  plane  wave. 
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While  these  assumptions  simplify  the  interaction  of  light  and  atmospheric 
particles,  they  do  not  always  represent  reality.  For  example,  many  particles  such  as  cubic 
salt  crystals  are  not  spherically  shaped.  However,  Mie  scattering  theory  can  still  be 
applied  in  that  case  since  the  irregular  particles  are  suspended  in  a  turbulent  media  and 
are  not  preferentially  oriented  in  one  direction  [21].  More  sophisticated  codes  extend  Mie 
scattering  to  include  spheres  with  layers  of  thin  films  for  various  purposes,  such  as 
modeling  water  layers  on  aerosol  particles  [22]. 

2.  Amplitude  Scattering  Matrix 

In  calculating  the  atmospheric  extinction  coefficients,  the  initial  goal  of  Mie 
scattering  theory  is  to  determine  the  scattered  field  amplitude  at  an  angle  0  and  a 
specified  distance  r  from  a  scattering  object,  in  this  case  a  suspended  particle  or  aerosol 
[9].  While  Mie  theory  utilizes  Maxwell’s  equations  to  determine  the  solutions  to  the 
fields,  both  inside  and  outside  of  the  scattering  object  or  scatterer,  HEL  propagation  only 
concerns  “far-field”  scattering,  where  r  »  /  /I  and  a  is  the  particle  radius.  In  this  case, 

the  fields  take  the  form  of  [9] 
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where  and  Ef  are  the  parallel  and  perpendicular  components  of  the  incident  electric 

field  and  and  are  the  respective  components  of  the  scattered  electric  field.  The 
wavenumber  k  —  2nlX  is  representative  of  light  in  air  since  the  refractive  index  of  air 


h  —  \.  The  amplitude  scattering  matrix 
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determines  the  amount  and 


polarization  of  scattered  light  at  a  given  angle  6  measured  with  respect  to  the  incident 
light. 


The  scattering  amplitudes  5*;  and  take  the  form  of  [23] 

“^1  (^)  =  (cos^)  (cos^)] , 
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where  n^  and  depend  only  on  the  scattering  angle  6  and  are  defined  as  [23] 
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In  Equations  (31)  and  (32),  P]  is  the  associated  Legendre  polynomial.  The  Mie 


coefficients  and  are  dependent  only  on  the  optical  wavelength  T,  the  scatterer 


radius  a  ,  and  the  scatterer  index  of  refraction  n  ,  and  are  defined  as  [23] 
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where  ^ and  denote  the  Hankel  functions  [19]. 


3.  Total  Cross  Sections 

The  total  scattering  cross  section  integrated  over  all  angles  takes  the  form  of  [9] 

n  2  cc 

'^.=5^Z(2<  +  1)(KI  +hr)-  (35) 

The  extinction  (scattering  +  absorption)  cross  section  is  given  by 

o  2  CO 

<7,=— X(2(  +  l)SH(a,+i,),  (36) 

LTl  n-i 

where  9i(. . .)  denotes  the  real  part.  Finally,  the  total  absorption  cross  section  is  [9] 

(37) 

which  is  simply  the  difference  between  the  extinction  and  scattering  cross  sections.  The 
absorption  and  extinction  coefficients  are  calculated  from  their  respective  cross-sections 
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by  integrating  over  the  contributions  from  different  particle  sizes.  For  example,  the 
aerosol  scattering  coefficient  [5^  is  calculated  by  integrating  [9] 

Pa{^)  =  \lN{r)c7^{r)dr ,  (38) 

where  Ni^r^dr  is  the  number  of  particles  with  a  radius  between  r  and  r  +  dr  per  unit 
volume. 


4.  Rayleigh  Scattering  Limit  for  Determining 


Molecular  scattering  assumes  that  particles  of  size  a  are  much  smaller  than  the 
wavelength  of  the  laser  light  (a  <s;  /i.) .  In  the  Rayleigh  scattering  limit,  the  scattering 

cross-section  depends  on  the  wavelength  of  electromagnetic  energy  as  qc  ,  so  shorter 
wavelengths  scatter  more  readily.  A  famous  example  of  this  principle  is  the  color  of  the 
daytime  sky:  because  blue  light  is  of  a  shorter  wavelength  than  red  light,  it  Rayleigh 
scatters  more  readily  and  thus  the  sky  appears  blue. 


In  this  case,  the  Rayleigh  scattering  approximation  can  be  utilized  to  determine 
the  scattering  cross  section  by  [9] 
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n  =\n\  . 

The  molecular  scattering  coefficient  can  then  be  calculated  from  the  cross  section 
following  the  same  method  as  Equation  (38). 


Figure  11  visualizes  the  general  scattering  cross  section-dependent  trends 
associated  with  Mie  scattering  theory.  In  the  low-frequency  regime  the  circumference  of 
an  arbitrary  spherical  particle  is  smaller  than  the  wavelength  of  energy  passing  over  it, 
and  thus  Rayleigh  scattering  (  oc  )  describes  the  scattering  cross  section  cr^ . 
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Figure  1 1 .  Mie  theory  scattering  behavior  with  respect  to  the  cross-section  of 
the  scatterer.  The  Rayleigh  scattering  limit  is  indicated  in  the  figure. 

5,  Aerosol  Absorption  and  Scattering  Coefficients  and 

Aerosol  sizes  are  typically  on  the  order  of  an  optical  wavelength  (a  « A) . 
According  to  Mie  theory,  aerosol-sized  particles  scatter  and  absorb  electromagnetic 
energy  in  a  manner  weakly  dependent  on  the  wavelength  of  light  passing  through  them 
[9].  Figure  11  visualizes  this  trend;  In  the  high  frequency  (optical)  limit  characterized  by 
the  resonance  when  a~  A,  the  cross  section  is  roughly  independent  of  the  wavelength. 
As  previously  stated,  the  aerosol  absorption  and  scattering  coefficients  are  calculated 
following  the  method  in  Equation  (38).  Flowever,  unlike  the  Rayleigh  scattering 
approximation  for  calculating  aerosol  coefficients  have  no  similar  approximating 
solutions. 

C.  MODTRAN 

There  are  several  well-established  radiative  transfer  codes  that  are  used  to 
generate  extinction  coefficients  needed  to  model  laser  propagation.  One  of  the  most 
widely  used  of  these  is  MODTRAN  (MODerate  resolution  atmospheric  TRANsmission) 
[25],  a  band-transmission  code,  which  utilizes  a  database  of  pre-calculated  coefficients 
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for  a  wide  variety  of  atmospherie  eonditions  and  aerosol  eontent.  MODTRAN  uses  the 
HITRAN  database  for  moleeular  absorption  speetra  [26]. 

MODTRAN  is  a  relatively  fast  (but  potentially  less  aeeurate)  method  of  obtaining 
optieal  properties  eompared  to  other  mueh  more  eomputationally  intensive  eodes  sueh  as 
the  line-by-line  LBLRTM  program  [27].  The  reason  line-by-line  eodes  are  more  CPU 
intensive  is  that  they  first  sum  over  all  absorption  lines,  then  add  eontinuum  absorption 
for  a  speeifie  wavelength,  then  ealeulate  Mie  absorption  and  seattering  eoeffioients  from 
physieal  prineiples  to  finally  determine  the  extinetion  eoeffieient. 

However,  while  MODTRAN  may  be  a  potentially  less  aeeurate  means  of 
generating  optieal  properties,  the  program  is  very  flexible  as  it  allows  for  various 
atmospherie  parameters  to  be  set  by  the  user. 

1,  Preset  Atmospheric  Models 

Table  2  lists  the  six  atmospherie  model  presets  that  MODTRAN  ean  use  to 
generate  atmospherie  optieal  properties.  The  user  ean  utilize  one  of  these  presets  as  is, 
modify  aspeets  of  these  presets,  or  speeify  user-defined  models. 


Atmospheric  Model  # 

Description 

1 

Tropical  Atmosphere  (15°  Latitude) 

2 

Mid-Latitude  Summer  (45°  North  Latitude) 

3 

Mid-Latitude  Winter  (45°  North  Latitude) 

4 

Sub-Arctic  Summer  (60°  North  Latitude) 

5 

Sub- Arctic  Winter  (60°  North  Latitude) 

6 

1976  US  Standard  Atmosphere 

Table  2.  MODTRAN  preset  atmospherie  models,  after  [25]. 
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The  six  preset  models  vary  in  several  eharaeteristies,  most  importantly  in  their 
temperature,  HjO,  and  O3  profdes.  Figure  12  displays  the  temperature  profdes  as 
funetions  of  altitude  for  eaeh  of  these  presets  [25]. 


Figure  12.  MODTRAN  model  atmosphere  temperature  profiles.  The  varying 
eolors  indieate  temperature  profiles  for  the  different  model 
atmospheres,  from  [25]. 


Naturally,  the  Sub-Aretie  Winter  preset  exhibits  the  eoldest  surfaee  temperature 
while  the  Tropieal  model  exhibits  the  warmest. 

Figure  13  displays  the  HjO  and  O3  profiles  as  funetions  of  altitude  for  the  six 
preset  atmospheres  [25]. 
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Figure  13.  MODTRAN  model  atmosphere  H2O  and  O3  profiles.  The  gaps  in  the 
data  are  the  result  of  the  breaks  in  the  vertieal  and  horizontal  axes,  as 

indieated,  from  [25]. 
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As  expeeted,  the  Tropieal  model  atmosphere  eontains  the  highest  eoneentration  of  H2O 
at  the  surfaee  and  the  Sub-Aretie  Winter  eontains  the  lowest.  The  O3  profde  is  slightly 
more  eomplieated,  with  all  models  eontaining  a  roughly  equal  eoneentration  at  ground 
level,  and  varying  greatly  at  the  peak  at  ~  20  km,  whieh  represents  the  ozone  layer. 


2,  Meteorological  Range  Presets 


MODTRAN  allows  for  varying  aerosol  densities  by  speeifying  the  ground 
meteorologieal  range  presets  at  23,  15,  10  or  5  km.  Meteorologieal  range  (MR),  in  this 
ease,  is  defined  by  Equation  (49)  as  [25] 

MR[km]  = - (41) 

EXT550[km-']  +  0.01159km-‘ 

where  EXT550[km'']  represents  the  surfaee  aerosol  extinetion  at  550  nm  and 


0.01159km"‘  is  the  surfaee  Rayleigh  seattering  eoeffieient  at  550  nm.  Meteorologieal 
range  ean  also  be  equated  to  the  more  subjeetive  visibility,  defined  as  the  maximum 
distanee  at  whieh  an  individual  ean  spot  and  identify  a  prominent  objeet  along  the 
horizon  at  day  (or  a  light  at  night),  by  the  approximate  Equation  (42) 


Visibility  « 


MR 
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(42) 
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If  no  meteorological  range  is  specified,  a  default  value  is  generated  from  presets 
associated  with  the  boundary  layer  chosen  in  the  particular  aerosol  model  (see  next 
section). 


3,  Aerosol  Models 

Another  option  within  MODTRAN  specifies  the  aerosol  model  preset  used  for  the 
surface  layer  (from  zero  to  two  kilometers  above  sea  level)  that  also  includes  a  default 
surface-meteorological  range.  While  MODTRAN  also  allows  the  user  to  specify  aerosol 
models  at  altitudes  above  the  surface  layer,  this  thesis  will  only  consider  engagements 
below  two  kilometers.  The  surface  layer  atmospheric  extinction  coefficient  dependence 
on  the  relative  humidity  is  determined  by  the  water-vapor  content  of  the  selected 
atmosphere  model.  Table  3  lists  the  primary  aerosol  presets  available  in  MODTRAN, 
with  model  #7  as  the  user-specified  aerosol  extinction  coefficient  option  excluded. 


Aerosol 

Model  # 

Description 

0 

No  aerosol  or  cloud  attenuation 

1 

RURAL  model;  Default  meteorological  range  set  to  23  km 

2 

RURAL  model;  Default  meteorological  range  set  to  5  km 

3 

NAVY  MARITIME  model;  Default  meteorological  range  set  based  on  wind  speed  &  relative 

humidity 

4 

MARITIME  model;  Default  meteorological  range  set  to  23  km 

5 

URBAN  model;  Default  meteorological  range  set  to  5  km 

6 

TROPOSPHERIC  model;  Default  meteorological  range  set  to  50  km 

8 

FOGl  model  (Advective  Fog);  Default  meteorological  range  set  to  0.2  km 

9 

FOG2  model  (Radiative  Fog);  Default  meteorological  range  set  to  0.5  km 

10 

DESERT  model;  Default  meteorological  range  set  by  wind  speed 

Table  3.  MODTRAN  main  aerosol  and  cloud  options,  after  [25]. 
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The  preset  most  commonly  used  in  this  study  was  the  #4  Maritime  model,  with  a 
default  ground  meteorological  range  of  23  km.  The  primary  purpose  of  the  Maritime 
model  is  to  study  the  propagation  of  infrared  energy  through  the  atmosphere  in  the  near 
sea  surface  layer,  an  area  in  which  all  Naval  vessels  would  operate  in,  making  the  model 
useful  for  this  research. 

This  model  utilizes  an  aerosol  profile  composed  of  a  sea-salt  component  and  a 
continental  component,  the  latter  of  which  is  primarily  comprised  of  water  soluble 
material  such  as  ammonium,  calcium  sulphate,  and  other  organic  compounds.  This 
aerosol  profile  does  not  contain  any  larger  dust  particles  due  to  the  assumption  that  these 
particles  would  rapidly  fall  out  of  the  atmosphere  as  a  result  of  gravitational  settling, 
which  occurs  on  the  order  of  hours  to  days  as  the  air  moves  from  land  and  across  the 
ocean  [21]. 

Furthermore,  the  Maritime  model  also  neglects  the  effects  of  fresh  sea-spray, 
which  occurs  in  the  lower  altitudes  above  the  sea  surface  at  10-20  meters  and  is  strongly 
dependent  on  surface  wind  speed.  Other  Maritime  aerosol  models  are  available  in 
MODTRAN  that  include  these  sea-spray  effects  and  could  represent  a  possible  extension 
to  this  research. 

4.  Cloud  and  Rain  Models 

Additional  options  specify  the  cloud  and  rain  presets  utilized  by  MODTRAN  in 
generating  atmospheric  extinction  coefficients.  The  rain  profiles  decrease  linearly  from 
the  ground  to  the  top  of  the  associated  cloud  model,  with  the  rain  cutting  off  at  the  cloud 
top  [20].  Table  4  lists  the  major  MODTRAN  cloud  and  rain  models,  with  preset  #  11  as 
the  user-specified  cloud  extinction  and  absorption  coefficient  option  excluded. 
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Rain  Model  # 

Description 

0 

No  clouds  or  rain 

1 

Cumulus  cloud  layer  with  base  at  0.66  km  and  top  at  3.0  km 

2 

Altostratus  cloud  layer  with  base  at  2.4  km  and  top  at  3.0  km 

3 

Stratus  cloud  layer  with  base  at  0.33  km  and  top  at  1 .0  km 

4 

Stratus/stratocumulus  cloud  layer  with  base  at  0.66  km  and  top  at  2.0  km 

5 

Nimbostratus  cloud  layer  with  base  at  0.16  km  and  top  at  0.66  km 

6 

2.0  mm/hr  ground  Drizzle  (Includes  ICLD  #  3  with  0.86  mm/hr  rain  at  1.0  km) 

7 

5.0  mm/hr  Light  rain  (Includes  ICLD  #  5  with  2.6  mm/hr  rain  at  0.66  km) 

8 

12.5  mm/hr  ground  Moderate  rain  (Includes  ICLD  #  5  with  6.0  mm/hr  rain  at  0.66  km) 

9 

25.0  mm/hr  ground  Heavy  rain  (Includes  ICLD  #  1  with  0.2  mm/hr  rain  at  3.0  km) 

10 

75.0  mm/hr  ground  Extreme  rain  (Includes  ICLD  #  1  with  1 .0  mm/hr  rain  at  3.0  km) 

18 

Standard  Cirrus  model  with  64  pm  mode  radius  for  ice  particles 

19 

Sub-visual  Cirrus  model  with  4  pm  mode  radius  for  ice  particles 

Table  4.  MODTRAN  cloud  and  rain  models,  after  [25]. 


D.  LEEDR 

LEEDR  (Laser  Environmental  Effects  Definition  and  Reference)  is  a  MATLAB- 
based  atmospheric  modeling  program  created  by  the  Air  Force  Institute  of  Technology 
(AFIT),  and  serves  a  similar  function  as  MODTRAN  in  characterizing  a  location’s 
atmosphere  based  on  various  user  inputs  such  as  location,  weather,  etc.  [28].  Because 
using  MODTRAN  to  generate  optical  parameters  for  model  atmospheres  containing 
variable  amounts  of  precipitation  proved  unsuccessful,  LEEDR  was  thus  implemented  to 
obtain  this  data  set.  All  LEEDR-generated  rainy  model  atmospheres  were  calculated 
using  atmosphere  types  and  aerosol  models  corresponding  to  those  available  in 
MODTRAN.  While  there  were  slight  variations  in  the  final  propagation  results  obtained 
from  the  model  atmospheres  generated  by  MODTRAN  and  LEEDR,  the  data  strongly 
agreed  at  the  1.064  pm  wavelength  (which  corresponds  to  the  nominal  laser  wavelength 
for  the  SSL). 
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V.  SIMULATION  METHODS 


In  order  to  accurately  and  efficiently  model  HEL  propagation  in  various 
atmospheric  conditions,  NFS  has  developed  ANCHOR  (Atmospheric  NFS  Code  for  HEL 
Optical  pRopagation).  ANCHOR  is  a  succinct  Matlab  script  that  utilizes  pre-calculated 
data  sets  containing  extinction  coefficient  values  across  a  specified  wavelength  spectrum. 
The  components  of  the  extinction  coefficient  -  namely  the  molecular  and  aerosol 
absorption  and  scattering  coefficients  (see  Chapter  III,  Section  A)  -  and  specified  laser 
parameters  are  subsequently  integrated  into  the  scaling  equations  governing  laser 
propagation.  The  on-target  laser  irradiance  and  power-in-the-bucket  values  are  obtained 
across  a  specified  wavelength  spectrum. 


Eour  specific  effects  govern  the  behavior  of  the  beam  as  it  propagates  to  the 
target:  Diffraction,  turbulence,  platform  jitter,  and  thermal  blooming.  The  linear  effects  of 
diffraction,  turbulence  and  jitter  are  assumed  to  produce  a  time-averaged  Gaussian 
irradiance  profile,  and  can  be  evaluated  independently  of  each  other.  These  effects  are  all 
implemented  within  the  “master  Equation  (43)  [29],”  which  describes  the  time-averaged 
irradiance,  or  intensity,  of  the  laser  beam  on  a  target  (in  W/m^). 


nkw. 


2~^  * 


(43) 


In  Equation  (43),  represents  the  total  output  power  from  the  laser  beam 


director,  is  the  time-averaged  radius  (1/e  in  the  intensity)  of  the  laser  spot  on 

target.  The  factor  comes  from  Beer’s  law  (see  Chapter  III,  Section  A),  where  e  is  the 

total  extinction  coefficient  due  to  atmospheric  absorption  and  scattering  (in  m"')  and  i  is 
the  distance  (in  m)  from  the  beam  director  to  the  target.  The  Sjg  factor  represents  the 

thermal  blooming  Strehl  ratio,  which  accounts  for  the  non-linear  reduction  in  the  on-axis 
irradiance  due  to  thermal  blooming. 


The  total  laser  beam  spot  size  on  target  is  determined  by  [29] 

2  2  2  2 

w  =  W ,  +  W  +  W . 


(44) 
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in  which  represents  the  root-mean-square  (rms)  spot  size  due  to  only  diffraction  in  a 

vacuum,  w,  the  rms  spot  size  due  to  only  turbulenee,  and  w-  the  rms  spot  size  due  to 

only  platform  jitter.  Sinee  these  effeets  are  independently  assumed  to  result  in  time 
averaged  Gaussian  irradianee  patterns,  their  eombined  rms  eontributions  add  in 
quadrature. 

A.  LINEAR  EFFECTS 

Before  the  linear  effeets  of  diffraetion,  turbulenee  and  jitter  governing  the 
behavior  of  a  propagating  laser  beam  can  be  implemented  into  Equation  (43),  they  have 
to  be  first  considered  and  quantified  independently. 


1. 


Vacuum  Diffraction: 


The  effects  of  laser  beam  diffraetion  in  a  vacuum  can  be  quantified  by  estimating 
the  beam’s  spot  size  due  to  diffraetion,  whieh  follows  as 

lU 


w,  ~  M- 
“  ttD 


(45) 


where  D  represents  the  beam  diameter  at  the  souree  and  M  represents  the  mode  faetor 
measuring  beam  quality.  If  the  laser  has  a  uniform  transverse  beam  shape,  then  D  is  also 
the  diameter  of  the  beam  direetor;  for  a  Gaussian  shape,  it  represents  the  diameter  of  the 
beam  at  the  1  /  e  point  in  the  irradianee.  The  mode  factor  is  a  unit-less  parameter  whieh 
represents  the  mode  eontent  of  a  multimode  laser  beam  obtained  from  a  stable  resonator. 
The  square  of  M ,  the  parameter,  eompares  the  produet  of  beam  divergenee  0  (in 
rad)  and  spot  size  tv  of  a  real  laser  beam  to  that  of  a  fundamental  Gaussian  referenee 
beam,  so  that  [30] 
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0  .w  , 

rej  rej 


(46) 


where  0^^^  represents  the  referenee  beam  divergence  and  the  reference  beam  spot 

size.  The  parameter  is  a  widely  used  quality  metric  for  laser  beams,  and  as  sueh  is 
often  defined  as  the  “beam  quality”  faetor.  However,  a  more  aeeurate  definition  of  is 
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the  “beam  propagation”  factor,  since  the  parameter  ultimately  provides  a  measure  of  how 
the  beam  will  propagate  through  free  space  [31].  For  an  ideal  Gaussian  beam,  =1, 

while  for  a  non-ideal  beam,  >  1 .  An  value  less  than  unity  is  not  possible. 


The  spot  size  due  to  diffraction  can  be  further  described  by  taking  into  account  the 
laser  beam’s  focal  range  F  ,  which  yields 

y  A 


2  2 


Fw 


(47) 


where  the  wavenumber  k  -Itt  /  A  and  represents  the  laser  beam  radius  at  the  source. 


For  this  research,  the  focal  point  of  the  laser  beam  was  placed  on  the  target,  so  F  =  £  and 
the  second  term  in  (47)  is  eliminated. 


2.  Atmospheric  Turbulence:  w, 


Turbulence  is  described  comprehensively  in  Chapter  III,  Section  B.  The  effect  of 
turbulence  on  the  laser  spot  size  on  target  can  be  estimated  by 
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(48) 


where  £  once  again  represents  the  distance  to  target  and  is  the  Fried  parameter,  which 


defines  a  length  scale  within  which  the  laser  beam  maintains  transverse  coherence  during 
its  propagation  path.  The  Fried  parameter  is  estimated  by 

/  T  \  -0.6 

r„«3.03 

V 


[— 1 


(49) 


where  characterizes  the  strength  of  the  turbulence,  as  discussed  in  Chapter  III, 
Section  B. 


3.  Platform  Jitter:  Wj 

Platform  jitter  is  defined  as  the  randomized  movement  of  the  laser  spot  position 
on  the  target  due  to  the  motion  of  the  laser  platform  as  well  as  tracking  error  by  the  active 
or  passive  tracking  system.  The  time-averaged  effect  of  jitter  can  be  quantified  by 
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considering  Equation  (50),  where  6^^^  defines  the  total  rms  angular  displacement  (in  rad) 
due  to  all  eontributions  from  jitter  and  t  is  the  distanee  (in  m)  to  target; 


w. 


e  L 


(50) 


B.  THERMAL  BLOOMING 


As  previously  discussed,  thermal  blooming  deseribes  the  non-linear  phenomenon 
in  which  ambient  air  surrounding  the  propagating  laser  beam  is  heated  by  the  beam  itself, 
resulting  in  a  ehange  in  the  loeal  air  density.  This  ehange  in  air  density  results  in  a  ehange 
in  the  refraetive  index  of  the  loeal  atmosphere,  ereating  a  “lens”  in  the  air  that  eauses  the 
beam  to  diverge.  The  direet  effeet  of  thermal  blooming  is  the  reduetion  of  the  laser  beam 
irradianee  on  target.  Thermal  blooming  is  diseussed  extensively  in  Chapter  III,  Seetion  C. 


1. 


Distortion  Number:  N , 


The  effeets  of  thermal  blooming  can  be  eharaeterized  by  the  thermal  distortion 
number  N ^  (Chapter  III,  Seetion  C.3),  whieh  is  defined  as  [4]; 
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(51) 


where  the  wavenumber  k-lTi/  A,  ~  1  kg  /  m^  represents  the  ambient  air  density, 

~10^  J/kg-K  is  the  speeifie  heat  of  air  at  eonstant  pressure  of  air,  <2(z)  is  the 

absorption  coeffieient  (in  m"'),  T(z)  =  represents  the  total  transmission  of  the  laser 

energy,  Hj.  (z)  =  dnl  dT  «  -10^®  K  '  is  the  change  in  the  refractive  index  with  respeet  to 

temperature,  v(z)  is  the  wind  speed  in  the  transverse  direetion  of  the  beam,  and  is 

the  beam  diameter.  In  all  of  the  above  expressions,  z  is  the  distance  from  the  beam 
direetor  and  is  measured  along  the  propagation  path. 


a.  Constant  Atmospheric  Properties 

Given  a  horizontal  beam  path  from  the  souree  to  the  target,  whieh  is  valid  only  if 
both  the  source  and  target  are  at  equal  altitudes  above  the  sea  level,  the  atmospheric 
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absorption  coefficient,  refractive  index,  and  wind  speed  are  all  assumed  constant  along 
the  propagation  path,  and  thus  the  distortion  number  N can  be  approximated  as 


b. 


Beam  Diameter:  D{^z) 


(52) 


The  beam  diameter  along  the  propagation  path  has  to  be  estimated  in  order 

to  evaluate  Equation  (51).  At  the  beam  director,  the  laser  beam  is  assumed  to  have  a 
diameter  of  .  At  the  target  (assumed  to  be  at  the  beam  waist),  the  beam  has  a  time- 

averaged  radius  defined  by  Equation  (44).  Assuming  that  the  beam  propagates  between 
these  two  points  and  the  fundamental  mode  is  roughly  Gaussian,  the  beam  diameter  along 
the  path  can  be  estimated  as 

+  ■  <53) 

In  (53),  represents  the  laser  beam  spot  radius  at  the  target,  z  is  the  distance  from  the 
target  and  is  the  “effective”  Rayleigh  range,  which  is  defined  such  that 

Z)(z  =  -f)  =  Z)o  [5],  and  visualized  in  Eigure  14.  Eollowing  this  definition,  the 
“effective”  Rayleigh  range  can  be  described  by 


if 


(54) 
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Figure  14. 


Visualization  of  laser  beam  profile  when  Z)(z 


Do  >  •  After  [6]. 


-1)  =  Dq  and 


However,  if  ,  then  the  waist  of  the  beam  is  loeated  at  the  beam  director 

(now  at  z  =  0)  and  the  target  is  located  at  z  =  f ,  which  is  visualized  by  Figure  (15).  The 
beam  diameter  is  then  calculated  as 

D{z)  =  D  +  ,  (55) 

V  / 

where 


if 


Do<^w,^f 


(56) 
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tot 


Figure  15.  Visualization  of  laser  beam  profile  when  target  is  located  at  z  =  f 

and  Dq  <  ,  after  [6]. 
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2. 


Strehl  Ratios:  S. 


The  thermal  blooming  Strehl  ratio  Sj.^  represents  a  ratio  of  the  peak  irradianee  of 
a  laser  beam  to  the  peak  irradianee  of  a  referenee  beam  where  no  thermal  blooming  has 
oeeurred.  This  ratio  is  utilized  by  the  ANCHOR  seript  to  ineorporate  thermal  blooming 
effeets  into  the  time-averaged  irradianee  and  power- in-the-bueket  ealeulations  [29].  The 
Strehl  ratio  for  a  Gaussian  laser  beam  shape  can  be  estimated  by  [4] 

Sg= - ^ - r,  (57) 

""  1  +  0.0625A^ 

while  the  Strehl  ratio  for  a  uniform,  or  “top-hat,”  laser  beam  shape  can  be  calculated  from 


^  1  +  O.OlA)," 

where  once  more  represents  the  thermal  distortion  number  defined  by  Equation  (51). 
For  the  propagation  simulations  in  this  research,  the  beam  shape  described  by  (58)  was 
utilized. 


C.  IRRADIANCE  LINEWIDTH  CONVOLUTION 

While  Equation  (43)  adequately  calculates  time-averaged  irradianee  profiles  for 
various  laser  beam  configurations  over  a  broad  wavelength  spectrum,  the  calculation  by 
itself  does  not  take  into  account  the  finite  linewidth  of  a  laser  (—0.1 -1.0%  for  HEEs). 
This  characteristic  is  visualized  in  Figure  16,  in  which  the  horizontal  axis  represents  a 
small  portion  of  the  electromagnetic  spectrum  and  the  vertical  axis  defines  the  strength  of 
molecular  absorption  for  each  wavelength.  The  red  line  indicates  a  typical  1%  HEE 
linewidth.  As  evidenced  by  the  figure,  many  detailed  absorption  features  are  contained 
within  the  linewidth,  and  the  overall  on-target  irradianee  profile  of  the  beam  is  greatly 
dependent  on  how  these  features  affect  the  laser. 
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Wavelength  (pm) 

Figure  16.  Absorption  features  contained  within  a  finite  1%  laser  linewidth.  The 
red  line  indicates  the  spectral  range  of  the  HEL,  after  [26]. 


In  order  to  more  realistically  model  laser  irradiance  profiles,  the  variations  in  the 
irradiance  due  to  the  absorption  features  within  a  finite  laser  linewidth  must  be  averaged. 
Considering  a  distribution  of  calculated  irradiance  profiles  across  a  specified  frequency 
spectrum,  a  typical  linewidth  can  be  described  by  the  change  in  laser  power  as  a  function 
of  the  change  in  frequency  dP  /  df ,  in  which 

dP 


[—df^P 
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(59) 


Combining  Equations  (59)  and  (43),  the  general  formula  for  the  time-averaged  irradiance 
for  an  arbitrary  linewidth  becomes 


■dP!  df 


nw. 
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(60) 


To  simplify  the  above  equation,  the  assumption  is  made  that  the  laser  frequency  spectrum 
has  a  uniform  (“top-hat”)  distribution  centered  at  frequency  and  with  a  linewidth  A/ , 

as  shown  in  Figure  17.  This  simple  top-hat  approximation  can  be  improved  upon  in  the 
future,  but  for  now  adequately  captures  important  spectral  features. 
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Figure  17.  Uniform  distribution  of  a  laser  line  centered  at  frequency  /„  with 

linewidth  A/ . 


Following  Equation  (43),  the  time-averaged  and  linewidth-averaged  irradiance 


then  becomes 


(61) 


In  implementing  this  calculation,  which  amounts  to  a  “moving  average”  of  irradiance 
values  across  a  full  frequency  spectrum.  Equation  (61)  is  discretized  in  order  to  yield 


(62) 


where  Sf  is  the  frequency  difference  between  adjacent  sampling  points  and  Af  /  f 
represents  the  aforementioned  linewidth.  Typically  the  linewidth  is  ~  0.1 -1.0%  for 
HEEs. 

D,  ADDITIONAL  HEL  PERFORMANCE  METRICS 

While  irradiance  on  target  is  an  important  parameter  to  evaluate,  it  is  not  the  only 
way  to  estimate  laser  performance  and  damage. 
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1. 


Power-in-the-bucket:  Pg 


Power-in-the -bucket  Pg  is  defined  as  the  amount  of  power  from  the  HEL  that 
falls  within  a  pre-defined  circle  (the  “bucket”)  with  radius  on  the  target.  Thermal 
blooming  causes  the  laser  beam  to  smear  out  in  a  complex  manner,  so  an  effective  beam 
size  is  estimated  as 


'TB 


(63) 


w. 


tot,eff 


where  wt,  is  determined  from  linear  effects.  The  irradiance  can  then  be  written  as 
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(64) 


If  Tg  >  ,  then  the  power-in-the-bucket  is  approximated  by 

=  (65) 

which  is  to  say  that  all  the  power  reaching  the  target  plane  falls  within  the  bucket.  If 
Vg  <  ,  the  power-in-the-bucket  can  be  simply  approximated  by 

Pb  ~  {l)[7rrl) .  (66) 

2.  Dwell  Time: 

A  more  physical  parameter  that  can  be  used  to  evaluate  laser  beam  performance  is 
the  dwell  time  for  melting  a  through  a  certain  amount  of  target  material.  Given  an 

arbitrary  material  of  a  specified  thickness,  one  could  estimate  the  length  of  time 
necessary  for  a  laser  beam  to  remain  on  target  in  order  to  melt  that  given  volume  of  the 
material. 


a.  Conductive  and  Radiative  Power  Transfer/Loss 

As  a  laser  beam  dwells  on  a  target,  a  certain  amount  of  the  laser’s  energy  is  lost 
due  to  energy  conduction  away  from  the  target  into  the  surrounding  material.  This 
conductive  power  transfer  loss  P^  (in  W)  can  be  quantified  by  [6] 
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(67) 


_lnrdK{T^-T^) 

^  L 

in  which  r  represents  the  radius  of  the  on-target  laser  spot  (in  m),  d  is  the  thickness  of 
the  material,  k  is  the  thermal  eonduetivity  of  the  material,  represents  the  melting 
temperature  of  the  target,  Tq  represents  the  ambient  and  initial  temperature,  and  L 
represents  the  distanee  through  the  material  over  whieh  the  temperature  transitions  from 
to  Tq  .  This  distanee  depends  on  the  material  properties  and  the  dwell  time;  for  typieal 

metals  and  dwell  times  of  ~  10  s,  L~2  em  [6]. 


The  second  power  loss  meehanism  is  radiative  power  loss  ,  whieh  oeeurs  when 

the  heated  target  area  radiates  as  a  black  body.  The  radiative  power  loss  can  be  estimated 
by  [6] 

P,=crsA{T,:-T:),  (68) 

where  (J  is  the  Stefan-Boltzmann  eonstant  (-5.670x10^*  W  m"^  K  "* ),  £  is  the 
emissivity  of  the  material  (dimensionless),  and  the  radiating  surfaee  area  A  -  2nr^ . 


The  total  power  loss  P^  is  thus  defined  as  [6] 


(69) 


In  order  to  caleulate  the  total  power  loss  ,  the  spot  size  r  must  first  be  estimated.  From 
Equation  (63),  the  effeetive  laser  spot  size  is  estimated  as 


r  =  w. 


tot,eff 


\J^TB 


(70) 


This  radius  depends  on  the  wavelength  of  the  laser,  atmospherie  eonditions  and  the 
platform/target  geometry. 


b.  Time  Required  to  Melt  the  Target 

Utilizing  the  previous  eomputations,  the  usable  melting  power  delivered  to  the 
target  ean  be  quantified  by 

P„u=P„-Pl.  (71) 
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in  which  represents  the  laser  power  delivered  to  the  target.  Finally,  if  >  0 ,  the 


dwell  time  can  be  quantified  by  [6] 


(72) 


in  whieh  represents  the  energy  required  to  melt  a  given  amount  of  the  target. 


The  dwell  time  depends  on  the  total  range  to  the  target  since  the  melting  power 
delivered  to  the  target  decreases  as  the  range  to  the  target  increases.  One  can  further  infer 
that  if  the  melting  power  delivered  to  the  target  does  not  overcome  the  power  loss  due  to 
conduction,  the  dwell  time  effectively  approaches  infinity,  in  which  case  the  laser  beam 
will  never  significantly  heat  the  target  material  to  its  melting  point. 


3.  Discussion 

To  evaluate  the  effectiveness  of  laser  beam  damage  to  a  target,  neither  the  beam 
irradiance  nor  the  power-in-the-bucket  ean  be  independently  considered.  In  fact,  both  are 
helpful  to  fully  characterize  the  damage.  The  time-averaged  irradiance  provides  valuable 
information  because  it  quantifies  the  maximum  intensity  of  a  laser  beam  on  a  speeified 
target  and  is  irrespective  of  the  size  of  the  target.  However,  the  irradiance  does  not 
indicate  by  itself  how  much  power  is  being  delivered  to  the  target  area. 

Conversely,  the  power-in-the-bucket  data  provides  a  measure  of  how  mueh  power 
emitted  by  a  laser  beam  is  confined  to  a  particular  spot  on  a  target.  Beeause  a  substantial 
spot  size  on  target  is  required  in  order  to  inflict  significant  damage,  P^  provides  valuable 

insight  into  the  lethality  of  the  laser  beam.  However,  this  metric  depends  on  the  size  of 
the  bucket  selected,  under  consideration. 

Dwell  time,  which  depends  on  the  material  properties  of  the  target  and  the  melting 
volume  needed  for  a  kill,  can  provide  guidance  for  operational  procedures  in  various 
engagement  scenarios.  Together,  all  three  metrics  complement  each  other  and  provide  a 
well-rounded  characterization  of  laser  performance. 
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E.  MODEL  VALIDATION 

In  order  to  ensure  that  the  ANCHOR  code  produces  realistic  laser  propagation 
results,  irradiance  values  at  specific  wavelengths  were  compared  to  the  irradiance  results 
obtained  from  Wavetrain.  Wavetrain  is  a  propagation  code  produced  by  MZA  Associates 
Corporation  based  in  Albuquerque,  NM  and  Dayton,  OH,  and  differs  from  the  ANCHOR 
script  as  it  actually  simulates  laser  beam  propagation  by  solving  the  paraxial  wave 
equation  [32],  This  method  is  self-consistent  and  potentially  more  accurate  than  the 
scaling  method  used  in  ANCHOR,  but  also  several  orders  of  magnitude  slower. 

Wavetrain  allows  for  irradiance  calculations  to  be  carried  out  in  model 
atmospheres  with  thermal  blooming  effects  turned  on  or  off,  and  similar  data  can  be 
obtained  from  the  ANCHOR  script.  Irradiance  values  from  Wavetrain  were  determined 
by  averaging  its  output  over  many  time  steps  after  the  pattern  reached  steady-state,  which 
was  ~  0.1  s. 

Table  5  displays  ANCHOR  and  Wavetrain  time-averaged  irradiance  profiles  at 
various  atmospheric  turbulence  values  taking  into  account  the  effects  of  turbulence  and 
thermal  blooming  for  the  /I  =  1.064  pm  wavelength,  which  corresponds  to  the  standard 
wavelength  for  a  SSL. 
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Value 

Data  Type 

Wavetrain  (W/m^) 

ANCHOR 

(W/m^) 

0 

Without  TB 

3.6x10* 

3.6x10* 

With  TB 

1.7x10* 

1.4x10* 

10l6 

Without  TB 

2.5x10* 

1.9x10* 

With  TB 

1.3x10* 

8.3x10’ 

lO'® 

Without  TB 

6.9x10^ 

2.5x10’ 

With  TB 

4.1x10^ 

1.3x10’ 

lo'^ 

Without  TB 

4.1x10^ 

1.7x10® 

With  TB 

4.9x10® 

1.2x10® 

Table  5.  Irradiance  propagation  data  comparison  for  1.064  pm  laser 

beam. 


Table  6  displays  the  identical  parameters  to  Table  5,  but  for  the  A  =  3.82  pm 
wavelength. 


C’  Value 

Data  Type 

Wavetrain  (W/m^) 

ANCHOR 

(W/m’) 

0 

Without  TB 

3.5x10’ 

3.5x10’ 

With  TB 

1.9x10’ 

1.9x10’ 

JO-16 

Without  TB 

3.4x10’ 

3.4x10’ 

With  TB 

1.9x10’ 

1.8x10’ 

10*^ 

Without  TB 

2.6x10’ 

2.2x10’ 

With  TB 

1.5x10’ 

1.2x10’ 

lO'"' 

Without  TB 

6.0x10® 

3.2x10® 

With  TB 

5.8x10® 

2.2x10® 

Table  6.  Irradiance  propagation  data  comparison  for  3.82  pm  laser 

beam. 
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As  evidenced  by  the  compared  data  sets,  at  small  turbulence  levels  the  ANCHOR 
and  Wavetrain  results  exhibit  negligible  deviation.  As  turbulence  increases  to 
cl  «  10"'"' m"^'^ ,  ANCHOR  and  Wavetrain  differ  by  a  factor  of  ~  two,  with  the  irradiance 

predicted  by  ANCHOR  consistently  smaller  between  the  two  codes.  On  closer  inspection, 
it  appears  ANCHOR  is  overestimating  the  degradation  of  intensity  due  to  turbulence 
relative  to  Wavetrain.  Even  so,  despite  the  differences  that  arise  between  the  codes  under 
these  conditions,  both  methods  predict  laser  irradiance  values  well  within  an  order  of 
magnitude  of  each  other.  It  should  be  noted  that  Wavetrain  takes  ~  one  minute  of 
computation  time  per  run;  in  that  same  time,  ACHOR  completes  ~  10^  runs. 

A  trend  in  the  data  arises  from  the  observation  that  as  the  strength  of  turbulence 
increases,  the  strength  of  thermal  blooming  decreases.  This  is  due  to  turbulence  in  the 
atmosphere  causing  the  laser  beam  to  spread  out,  lowering  the  peak  irradiance  of  the 
beam  and  preventing  thermal  blooming  from  taking  effect. 
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VI.  RESULTS 


The  simulated  laser  propagation  results  obtained  from  ANCHOR  were  analyzed 
and  eatalogued  in  various  plots.  Due  to  the  effieiency  of  the  ANCHOR  seript,  each 
generated  data  plot  represents  tens  to  hundreds  of  thousands  of  simulation  runs;  nearly  10 
million  runs  in  total  were  conducted  for  this  thesis. 

A,  STANDARD  MODEL  PARAMETERIZATION 

To  effectively  utilize  the  ANCHOR  code  in  evaluating  laser  performance, 
standard  laser  and  atmosphere  models  were  chosen  in  order  to  establish  a  baseline 
performance  profile  for  the  various  studies  conducted.  Table  7  displays  the  laser  and 
target  parameters  that  define  the  standard  model  utilized  in  the  ANCHOR  simulations. 


Laser  Output  Power  P 

100  kW 

Laser  Wavelength  X 

0.5  pm  -  5  pm 

Laser  Linewidth  A/7  / 

1% 

On-Target  Bucket  Size 

5  cm 

Beam  Quality  M~ 

3 

Beam  Diameter  at  Source  D 

0.5  m 

Platform  Jitter  6 

rms 

2.0x10^  rad 

Target  Range 

5000  m 

Engagement  Geometry 

Horizontal;  10  m  Above  Water  Surface 

Adaptive  Optics 

None 

Table  7.  Standard  model  laser  and  target  parameters  for  ANCHOR. 
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Table  8  displays  the  atmospherie  parameters  that  define  the  standard  model 
utilized  in  the  ANCHOR  simulations. 


Atmosphere  Type 

Tropical,  15“  Latitude  [MODTRAN  Model  #1] 

Aerosol  Composition 

Maritime  [MODTRAN  Model  #4] 

Temperature,  Tq 

300  K  (80“  F) 

Absolute  Humidity 

19.7  g/m^  (80%  Relative  Humidity) 

Meteorological  Range 

23  km  (Clear  Day) 

Precipitation 

No  Rain 

Turbulence  C“ 

ft 

(Moderate) 

Wind 

Transverse,  10  m/s 

Table  8.  Standard  model  atmospherie  parameters  for  ANCHOR. 


Finally,  Figure  18  displays  a  simple  visualization  of  the  standard  model  target 
engagement  seenario,  in  whieh  both  the  laser  platform  and  target  are  positioned  10  m 
above  the  water  surface  and  are  separated  by  a  range  of  5000  m.  Once  again,  the  diameter 
of  the  laser  beam  at  the  source  (and  thus  the  beam  director  diameter)  is  0.5  m. 


Figure  18.  Visualization  of  standard  model  target  engagement  scenario. 
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B. 


STANDARD  MODEL  STUDIES  AT  5  KM  TARGET  RANGE 


Evaluations  of  the  standard  model  were  conducted  by  first  selecting  a  single 
parameter  (such  as  turbulence,  haze,  output  power,  etc.)  and  varying  it  over  a  realistic 
range  while  keeping  the  other  parameters  fixed  to  the  standard  model.  The  associated 
irradiance  and  power-in-the-bucket  profiles  were  then  determined  for  each  simulation. 
This  method  allowed  various  known  trends  associated  with  laser  propagation  to  be 
verified  and  visualized. 

In  order  to  ground  the  visualized  laser  performance  to  realistic  operational 
requirements,  a  “melt  threshold”  is  shown  in  the  irradiance  plots.  This  melt  threshold 
represents  the  minimum  on-target  laser  irradiance  required  to  melt  a  certain  volume  of 
the  target  material,  and  is  a  function  of  the  desired  area  to  damage  on  the  target  and  the 
laser  dwell  time.  In  this  case,  the  target  was  assumed  to  be  an  arbitrary  aircraft 
(unmanned  air  vehicle,  helicopter,  etc.)  with  a  body  comprised  of  1  mm- thick  aluminum. 
In  order  to  attain  the  melt  threshold,  the  laser  must  deliver  enough  energy  to  the  target  for 
the  aluminum  to  reach  its  melting  temperature  and  continue  to  deliver  enough  energy  to 
melt  through  the  skin,  all  within  a  10-second  laser  dwell  time. 


Table  9  displays  the  various  properties  of  aluminum  necessary  for  determining 
this  melt  threshold. 


Density  p 

2700  kg/m^ 

Melting  Temperature 

933  K 

Heat  Capacity 

897  J/kg-K 

Latent  Heat  of  Fusion  H 

400  kJ/kg 

Thermal  Conductivity  k 

237  W/m-K 

Emissivity  s 

0.05 

Table  9.  Physical  properties  of  aluminum,  from  [33]. 
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Estimating  the  melt  threshold  requires  knowledge  of  the  target  material  mass  m 
and  the  total  energy  neeessary  to  melt  this  mass  ,  whieh  is  eomprised  of  the  energy 

required  to  reaeh  the  material’s  melting  temperature  and  the  energy  required  to  melt 
the  target  onee  the  melting  temperature  is  reaehed  .  Assuming  the  desired  region  to 

melt  is  eireular  with  a  radius  r  -  0.05  m  and  material  thiekness  d  -  0.001m,  the  mass  m 
of  the  targeted  material  ean  be  estimated  by 

m  =  pnr^d  =  0.021  kg.  (73) 

Next,  the  energy  required  to  heat  the  material  is  [6] 

a-='”C,(7;.-r.)  =  12.0kJ,  (74) 

where  =  300  K  is  the  ambient  temperature.  The  energy  required  to  melt  the  aluminum 
onee  the  melting  temperature  has  been  reaehed  is 

a..=«^  =  8.4kJ.  (75) 

Finally,  the  total  energy  needed  to  melt  is  the  sum  of  the  two  energy  values, 

Q,=a-+a.,,=2')-4U.  (76) 

With  knowledge  of  the  total  energy  requirement,  the  on-target  laser  irradianee 
required  to  meet  the  melt  threshold  ean  be  estimated  by  [6] 

( Q  3 

—  +  (Pc+P,) 

K.„  =  — - ^ - -  =  C-Se  MW/m",  (77) 

nr 

where  is  the  laser  beam  dwell  time  (diseussed  in  Chapter  IV,  Seetion  D.2),  is  the 
eonduetive  power  loss  from  Equation  (67)  and  is  the  radiative  power  loss  from 

Equation  (68).  For  the  irradianee  plots,  the  melt  threshold  will  be  identified  by  a  solid, 
horizontal  red  line  above  whieh  the  laser  will  likely  infliet  signifieant  damage  to  the 
target. 

The  following  plots  present  either  irradianee  or  power-in-the-bueket  data  on  the 
vertieal  axis  for  a  wavelength  range  of  0. 5-5.0  pm,  with  the  varying  eolors  representing 
the  different  values  of  the  parameter  being  ehanged.  Unless  otherwise  stated,  the  laser 
and  atmosphere  parameters  eonform  to  the  values  listed  in  Tables  7  and  8. 
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1,  Turbulence  Study  with  a  5  km  Target  Range 

Figure  19  displays  plots  of  the  laser  irradiance  for  the  standard  model,  where  the 

2  16  2/3  1 3 

C„  parameter  is  varied  from  10'  m'  (weak  turbulenee,  indicated  by  dark  blue)  to  10' 
m'  (strong  turbulence,  indicated  by  light  blue).  As  evidenced  by  the  plot,  turbulence 
corresponding  to  Cl  >  lO'^'^m”^'"^  may  prevent  a  successful  engagement  at  a  5  km  range, 

since  the  irradiance  generally  falls  near,  or  below,  the  required  melt  threshold  to  inflict 
significant  damage  to  the  target.  Furthermore,  turbulence  tends  to  affect  the  shorter 
wavelengths  more  than  longer  wavelengths  (cf  Equation  4),  a  trend  reproduced  by 
ANCHOR  and  visible  in  this  figure. 


Standard  Model  Turbulence  Study 

melt  threshold  .  1e-13m''-2/3  1e-14m''-2/3  1e-15m^-2/3  16-16  m^>2/3 


1  2  3  4  5 

Wavelength  (pm) 


Figure  19.  Standard  model  turbulence  study:  Irradiance  results  indicate 
successful  engagement  at  turbulences  weaker  than  Cl  -  lO’^'^m”^^^ . 

Figure  20  displays  plots  of  the  power- in-the-bucket  for  a  bucket  radius  of  5  cm  for 
the  same  study  as  Figure  19.  Notice  that  moderate  turbulence  (C^  =10”'^m“^^^)does  not 

decrease  power-in-the-bucket,  while  stronger  turbulence  tends  to  increase  overall  on- 
target  laser  spot  size.  The  moderate  turbulence  here  does  not  increase  the  spot  size 
enough  to  where  it  falls  outside  the  bucket,  and  thus  all  of  the  laser  power  remains 
confined  within  it.  Once  the  turbulence  reaches  extreme  levels,  however,  the  on-target 
spot  size  becomes  larger  than  the  bucket  size,  and  thus  the  power-in-the-bucket  drops. 
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Figure  20.  Standard  model  turbulenee  study:  Power-in-the-bueket  results 
indieate  signifieant  loss  of  bueket  power  in  the  presenee  of  strong 

turbulenee. 

2,  Variation  of  the  Laser  Beam  Diameter  at  the  Source  with  a  5  km 
Target  Range 

The  following  study  examines  the  trends  in  on-target  laser  irradianee  and  power- 
in-the-bueket  assoeiated  with  varying  the  laser  beam  diameter  at  the  source  at  a  target 
range  of  5  km  and  in  the  presence  of  weak  and  strong  turbulence. 

a.  Weak  Turbulence 

Figure  21  displays  plots  of  the  laser  irradianee  for  the  standard  model,  where  the 
beam  diameter  at  the  source,  and  thus  the  beam  director  size  D ,  has  been  varied  from 
0.3  m  (light  blue)  to  0.7  m  (dark  blue)  within  a  weak  turbulence  (C^  =10“''^m”^^^) 

environment.  In  weak  turbulence,  increasing  the  beam  director  size  shows  marked 
improvement  in  the  target  irradianee.  Furthermore,  shorter  wavelengths  benefit  less  from 
the  increase  in  the  beam  director  size  since  turbulence  and  jitter  contributions  tend  to 
dominate  the  time-averaged  spot-size. 
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1  2  3  4  5 


Wavelength  (gm) 

Figure  21 .  Standard  model  weak  turbulence  beam  director  size  study:  Irradiance 
results  indicate  marked  improvement  in  laser  performance  by 
increasing  beam  director  size  in  a  weakly  turbulent  environment. 


Figure  22  represents  the  power-in- the-bucket  results  for  the  same  study  as  Figure 
21.  Increasing  the  beam  director  size  has  little  effect  on  the  power- in-the-bucket  in  the 
presence  of  weak  turbulence  for  shorter  wavelengths;  for  longer  wavelengths  (>2  pm), 
the  reduction  of  the  spot  size  due  to  the  larger  beam  director  allows  more  of  the  power  to 
be  contained  within  the  bucket. 


standard  Model  Source  Laser  Beam  Diameter  Study  -  Weak  T urbulence  -  PIB 


Figure  22.  Standard  model  weak  turbulence  beam  director  size  study:  Power-in- 
the-bucket  results  indicate  little  effect  of  beam  director  size  on  power 
delivered  to  bucket  in  weakly  turbulent  environment. 
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b.  Strong  Turbulence 

Figure  23  displays  plots  of  the  laser  irradiance  for  the  standard  model,  in 
which  the  laser  beam  diameter  at  the  source,  and  thus  the  beam  director  size  D ,  has  been 
varied  from  0.3  m  to  0.7  m  within  a  strong  turbulence  (C^  =10“''^m“^^^)  environment.  It 

appears  that  increasing  the  beam  director  size  in  this  case  has  little  effect  on  the 
laser  performance,  since  the  strong  turbulence  contribution  dominates  laser  irradiance 
regardless  of  the  size  of  the  laser  beam  at  the  source. 


standard  Model  Source  Laser  Beam  Diameter  Study  -  Strong  Turbulence 


Figure  23.  Standard  model  strong  turbulence  beam  director  size  study: 

Irradiance  results  indicate  little  effect  of  beam  director  size  on  laser 
performance  in  turbulent  environment. 


Figure  24  displays  plots  of  the  power-in-the-bucket  results  for  the  same  study  as 
Figure  23.  Increasing  the  beam  director  size  has  once  again  little  effect  on  the  power 
delivered  within  the  bucket  on  the  target  within  a  strong  turbulence  environment. 
However,  in  comparison  to  Figure  22,  the  bucket  power  has  markedly  dropped  for  all 
wavelength  windows;  this  occurs  because  in  this  more  turbulent  environment  the  spot 
size  at  the  target  is  now  larger  than  the  bucket  size  regardless  of  the  wavelength. 
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standard  Model  Source  Laser  Beam  Diameter  Study  -  Strong  Turbulence  -  PIB 


Figure  24.  Standard  model  strong  turbulence  beam  director  size  study:  Power- 
in-the-bucket  results  indicate  little  effect  from  beam  director  size  on 
power  delivered  to  a  bucket  in  a  turbulent  environment. 

3,  Meteorological  Range  (Haze)  Study  with  a  5  km  Target  Range 

Figure  25  displays  plots  of  the  laser  irradiance  for  the  standard  model,  where  the 
haze  (meteorological  range)  has  been  varied  from  5  km  visibility  (hazy  day,  indicated  by 
light  blue)  to  23  km  (clear  day,  indicated  by  dark  blue).  Strong  haze  (visibility  <  5  km) 
prevents  a  successful  engagement  at  a  5  km  range,  as  the  melt  threshold  is  not  likely  to  be 
met  at  such  limited  visibility.  Furthermore,  haze  tends  to  affect  shorter  wavelengths  more 
than  longer  wavelengths.  As  such,  wavelengths  greater  than  3  pm  experience  less  of  a 
reduction  in  on-target  irradiance  than  the  0.5  pm  to  2.5  pm  wavelength  range.  This  is  a 
known  trend  reproduced  by  these  results  [5]. 


standard  Model  Haze  (Meteorological  Range)  Study 

-  mell  threshold  5  km  10  km  ■■■  15  km  23  km 

10®  j- - ■ - ^ ^  ^ — >  — - -  — - - — • — ^ — • 


1  2  3  4  5 

Wavelength  (pm) 


Figure  25.  Standard  model  meteorological  range  study:  Irradiance  results 
indicate  successful  target  engagement  at  visibilities  greater  than  5  km. 
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Figure  26  displays  plots  of  the  power-in-the-bucket  data  for  the  same  study  as 
Figure  25.  From  the  data,  it  appears  that  the  reduction  in  power-in-the-bucket  is 
proportional  to  the  reduction  in  visibility.  Instead  of  affecting  the  laser  spot  size,  like 
turbulence  does,  haze  instead  attenuates  the  laser  beam  as  it  propagates  to  the  target  via 
atmospheric  absorption  and  scattering. 


Figure  26.  Standard  model  meteorological  range  study:  Power-in-the-bucket 
results  indicate  overall  reduction  in  power  due  to  atmospheric 

attenuation. 


4,  Precipitation  Study  with  a  5  km  Target  Range 

Figure  27  displays  plots  of  the  laser  irradiance  for  the  standard  model,  where 
various  amounts  of  precipitation  have  been  introduced.  Even  the  presence  of  light 
precipitation  of  5.0  mm/hr  rain  prevents  the  likelihood  of  a  successful  engagement  at  a  5 
km  range  due  to  the  attenuation  of  the  laser  power  by  the  rain  droplets.  Irradiance 
estimates  in  heavier  precipitation  rates,  while  successfully  calculated,  fall  below  the  scale 
of  this  plot  at  the  5  km  target  range.  However,  engagements  at  shorter  ranges  may  still  be 
viable,  as  will  be  seen  later. 
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standard  Model  Rain  Study  >  5  km  Target  Range 


melt  threshold 

25.0  mm/hr  Heavy  Rain 

12.5  mm/hr  Moderate  Rain  ' 

5.0  mnVhr  Light  Rain 

No  Rain 

Wavelength  (pm) 


Figure  27.  Standard  model  precipitation  study:  Irradiance  results  indicate  that 
even  light  rain  prevents  likelihood  of  a  successful  engagement  at  a  5  km 

target  range. 


Figure  28  displays  plots  of  the  power-in-the-bucket  data  for  the  same  study  as 
Figure  27.  Again,  the  results  for  heavier  rain  rates  fall  below  this  scale  shown  in  the  plot. 
A  successful  engagement  is  unlikely  for  even  the  slightest  amount  of  precipitation  at  a  5 
km  target  range. 
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Figure  28.  Standard  model  precipitation  study:  Power-in-the-bucket  results 

indicate  significant  power  reduction  in  the  presence  of  light  rain  at  5  km 

target  range. 
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5,  Beam  Quality  Study  with  a  5  km  Target  Range 

Figure  29  displays  plots  of  the  laser  irradiance  profiles  for  the  standard  model.  In 
this  figure  the  beam  quality  factor  has  been  varied  from  1  (perfect  beam  quality, 
indicated  by  dark  blue)  to  10  (poor  beam  quality,  indicated  by  light  blue).  Longer  laser 
wavelengths  experience  a  greater  reduction  in  performance  with  degraded  beam  quality. 
The  reason  is,  once  again,  that  longer  wavelength  spot  sizes  are  dominated  by  vacuum 
diffraction  (cf.  Equation  45),  where  the  beam  quality  plays  a  significant  role. 


Standard  Model  Beam  Quality  Study 


Wavelength  (pm) 


Figure  29.  Standard  model  beam  quality  study:  Irradiance  results  indicate  a 
greater  reduction  in  laser  performance  at  longer  wavelengths  with 
decreased  beam  quality. 


Figure  30  displays  plots  of  the  power-in- the-bucket  data  for  the  same  study  as 
Figure  29.  Similar  to  the  irradiance  results,  degradation  in  the  beam  quality  causes  a 
reduction  in  power-in-the-bucket  for  longer  laser  wavelengths,  again  due  to  diffraction. 
The  increase  in  spot  size  reaches  a  point  to  where  some  of  the  laser  power  falls  outside 
the  bucket,  correlating  to  the  decreased  power  delivered  to  the  desired  target  area. 
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standard  Model  Beam  Quality  Study  -  PIB 


Figure  30.  Standard  model  beam  quality  study:  Power-in-the-bucket  results 
indieate  greater  reduction  in  bucket  power  at  longer  laser  wavelengths 
for  decreased  beam  quality. 

6,  Platform  Jitter  Study  with  a  5  km  Target  Range 

Figure  31  displays  plots  of  the  laser  irradiance  profiles  for  the  standard  model,  in 
which  the  platform  jitter  0^^^  has  been  varied  from  0  prad  (no  jitter,  indicated  by  dark 
blue)  to  10  prad  (extreme  jitter,  indicated  by  light  blue).  An  increase  in  platform  jitter 
results  in  a  reduction  of  laser  performance  for  all  wavelengths;  however,  shorter 
wavelengths  tend  to  be  affected  more.  While  the  absolute  contribution  of  jitter  to  the  spot 
size  is  wavelength  independent  (cf.  Equation  50),  the  relative  contribution  is  greater  for 
shorter  wavelengths  (cf.  Equation  44). 


standard  Model  Platform  Jitter  Study 


Figure  3 1 .  Standard  model  platform  jitter  study:  Irradiance  results  indicate  an 
equal  reduction  in  laser  performance  with  increased  platform  jitter  for 
all  wavelengths,  as  expected. 
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Figure  32  displays  plots  of  the  power-in-the-bueket  data  for  the  same  study  as 
Figure  31.  Similar  to  Figure  31,  an  increase  in  platform  jitter  increases  the  on-target  laser 
spot  size,  which  causes  more  laser  power  to  fall  outside  the  desired  target  area  (the 
bucket),  and  thus  less  power  is  delivered  within  the  bucket. 
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Figure  32.  Standard  model  platform  jitter  study:  Power-in-the -bucket  results 
indicate  a  reduction  in  laser  performance  with  increased  platform  jitter 

for  all  wavelengths. 

7.  Output  Power  Study  with  a  5  km  Target  Range 

Figure  33  displays  plots  of  the  laser  irradiance  results  for  the  standard  model,  in 
which  the  laser  output  power  P  has  been  varied  from  30  kW  to  200  kW.  From  the 
results,  it  appears  that  increasing  laser  output  power  from  30  kW  to  100  kW  results  in  a 
marked  performance  increase  of  the  laser.  However,  doubling  the  output  power  from  100 
kW  to  200  kW  has  little  benefit,  as  the  increased  laser  power  also  increases  thermal 
blooming,  which  causes  a  further  reduction  of  on-axis  intensity.  Again,  this  is  a  known 
trend  that  is  reproduced  by  the  ANCHOR  code  [5]. 
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Figure  33.  Standard  model  output  power  study:  Irradianee  results  indieate  a 
marginal  inerease  in  laser  performanee  when  doubling  laser  output 
power  from  1 00  kW  to  200  kW. 


Figure  34  displays  plots  of  the  power-in-the-bucket  data  for  the  same  study  as 
Figure  33.  While  doubling  the  laser  output  power  from  100  kW  to  200  kW  eauses  the 
effeets  of  thermal  blooming  to  be  more  pronouneed  (lowering  the  on-target  laser 
intensity),  power-in-the-bueket  appears  to  be  less  affeeted  by  this  phenomenon.  While 
thermal  blooming  eauses  the  on-target  laser  spot  size  to  deform  mueh  like  the  depietion 
in  Figure  6,  the  spot  appears  to  remain  generally  eonfined  within  the  desired  on- target 
bucket,  effectively  delivering  more  power  to  the  target  at  higher  laser  output  powers  and 
especially  for  shorter  wavelengths. 


Figure  34.  Standard  model  output  power  study:  Power-in-the-bucket  results 
indicate  reasonable  increase  in  laser  performance  by  doubling  output 
power  from  1 00  kW  to  200  kW. 
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c. 


STANDARD  MODEL  STUDIES  AT  VARIABLE  TARGET  RANGES 


The  previous  studies  eoneerned  the  effeets  of  various  atmospherie  eonditions  and 
laser  properties  on  overall  performanee  while  engaging  a  target  at  a  fixed  5  km  range. 
However,  this  approaeh  is  not  wholly  realistie  as  an  adversary  will  rarely  remain 
stationary  in  an  operational  environment.  The  following  studies  attempt  to  develop 
understanding  of  how  atmospherie  eonditions  affeet  laser  lethality  as  the  target  eloses  in 
on  the  laser  souree. 

1.  Target  Range  Variation  Study 

Figure  35  displays  plots  of  the  laser  irradianee  results  for  the  standard  model,  in 
whieh  the  target  range  has  been  varied  from  1000  m  (indieated  by  dark  blue)  to  7000  m 
(indieated  by  light  blue).  It  appears  that  on  a  elear  day  with  moderate  turbulenee  and  a  10 
m/s  erosswind,  a  target  is  likely  to  be  sueeessfully  engaged  at  ranges  eloser  than  7000  m 
for  systems  operating  at  the  nominal  ~1  pm  SSL  wavelength. 


Figure  35.  Standard  model  target  range  study:  Irradianee  results  indicate 
likelihood  of  a  successful  engagement  at  ranges  closer  than  7  km  in 
clear  weather  conditions. 


Figure  36  displays  plots  of  the  power-in- the-bucket  data  for  the  same  study  as 
Figure  35.  Similar  to  the  conclusions  drawn  from  Figure  35,  as  a  target  approaches  the 
laser  source,  the  laser  lethality  is  increased.  However,  as  the  target  closes  in  from  7000  m 
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to  1000  m,  the  power- in-the -bucket  increases  since  less  of  the  beam  is  attenuated  by  the 
atmosphere.  At  a  range  of  -1000  m,  nearly  all  of  the  power  from  the  source  falls  within 
the  bucket,  since  the  spot  size  is  relatively  small  even  for  the  longer  wavelengths. 
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Figure  36.  Standard  model  target  range  study:  Power-in-the-bucket  results 
indicate  increase  in  laser  performance  as  target  closes  in  on  laser 

source. 

2,  Precipitation  Target  Range  Studies 

The  following  irradiance  and  power-in-the-bucket  studies  evaluate  laser 
performance  at  various  target  ranges  and  in  the  presence  of  precipitation  of  increasing 
severity. 

a.  Drizzle 

Figure  37  displays  plots  of  the  laser  irradiance  results  for  the  standard  model  in 
the  presence  of  drizzle  (2.5  mm/hr  rain),  with  a  variation  in  target  proximity  from  1000  m 
to  5000  m.  While  laser  lethality  in  the  presence  of  drizzle  is  marginal  at  best  for  a  target 
range  of  5  km,  the  on-target  irradiance  jumps  by  more  than  an  order  of  magnitude  as  the 
target  closes  in  to  3000  m,  effectively  making  the  weapon  highly  lethal  against  typical 
targets  at  this  range. 
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Figure  37.  Standard  model  drizzle  target  range  study:  Irradiance  results  indicate 
likelihood  for  successful  engagement  at  ranges  closer  than  5  km. 

Figure  38  displays  plots  of  power- in-the -bucket  results  for  the  same  study  as 
Figure  37.  Similar  to  the  results  from  Figure  37,  as  a  target  approaches  the  laser  source, 
the  power  delivered  to  the  bucket  on  the  adversary  increases  for  ah  laser  wavelengths. 
This  increase  in  power-in-the -bucket,  however,  stems  from  the  fact  that  the  atmospheric 
precipitation  attenuates  laser  power  very  effectively.  The  target  has  to  be  in  closer 
proximity  to  the  source  to  be  hit  by  a  substantial  amount  of  laser  power. 


standard  Model  -  Drizzle  Rain  (2.5  mm/hr)  Study  -  PIB 
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Figure  38.  Standard  model  drizzle  target  range  study:  Power- in-the -bucket 
results  indicate  marked  improvement  in  laser  performance  at  target 
ranges  closer  than  5  km. 
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b.  Light  Rain 

Figure  39  displays  plots  of  the  laser  irradiance  results  for  the  standard  model  in 
the  presence  of  light  rain  (5.0  mm/hr  rain),  with  a  variation  in  target  proximity  from  1000 
m  to  5000  m.  Similar  to  the  results  for  2.5  mm/hr  rain,  stronger  precipitation  reduces  the 
on-target  laser  intensity  by  attenuating  the  laser  power.  The  target  needs  to  be  within 
3000  m  for  the  laser  to  effectively  engage  it. 


Figure  39.  Standard  model  light  rain  target  range  study:  Irradiance  results 
indicate  likelihood  of  successful  target  engagement  at  ranges  closer 

than  3  km. 


Figure  40  displays  plots  of  the  power-in-the-bucket  data  for  the  same  study  as 
Figure  39.  Similar  to  the  power-in-the-bucket  results  for  2.5  mm/hr  drizzle,  a  target  range 
of  less  than  2000  m  to  the  source  results  in  an  increase  of  power-in-the-bucket  by  almost 
an  order  of  magnitude. 
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Figure  40.  Standard  model  light  rain  target  range  study:  Power-in-the-bucket 
results  indicate  marked  increase  in  laser  performance  for  target  ranges 

closer  than  3  km. 


c.  Moderate  Precipitation 

Figure  41  displays  plots  of  the  laser  irradiance  results  for  the  standard  model  in 
the  presence  of  moderate  precipitation  (12.5  mm/hr  rain),  in  which  the  range  to  target  has 
been  varied  from  1000  m  to  4000  m.  For  shorter  wavelengths  the  laser  irradiance  exceeds 
the  melt  threshold  at  a  target  range  of  2  km.  At  a  1  km  range,  the  on-target  laser  intensity 
exceeds  the  melt  threshold  by  a  large  factor. 


Figure  41 .  Standard  model  moderate  rain  target  range  study:  Irradiance  results 
indicate  likelihood  of  a  successful  engagement  within  a  target  range  of 

2  km. 
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Figure  42  displays  plots  of  the  power-in-the-bueket  data  for  the  same  study  as  Figure 
41.  Similar  to  previous  preeipitation  results,  the  power-in-the-bucket  inereases 
proportionally,  and  relatively  equally  for  all  wavelengths,  to  the  target  proximity  from  the 
source.  However,  the  laser  only  truly  becomes  effective  for  target  ranges  closer  than  2000  m. 


Standard  Model  -  Moderate  Rain  (12.5  mm/hr)  Target  Range  Study  -  PIB 
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Figure  42.  Standard  model  moderate  rain  target  range  study:  Power-in-the- 
bucket  results  indicate  marked  increase  in  laser  performance  for  target 

ranges  closer  than  2  km. 


d.  Heavy  Precipitation 

Figure  43  displays  plots  of  the  laser  irradiance  results  for  the  standard  model  in 
the  presence  of  heavy  precipitation  (25.0  mm/hr  rain),  and  as  the  target  range  is  varied 
from  1000  m  to  3000  m.  The  on-target  laser  irradiance  exceeds  the  melt  threshold  at  a  1 
km  target  range.  Thus,  even  in  heavy  rain,  a  100  kW-class  HEL  is  a  viable  weapon 
against  targets  within  1000  m  of  the  source. 
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Figure  43.  Standard  model  heavy  rain  target  range  study:  Irradianee  results  indieate 
likelihood  of  a  sueeessful  engagement  at  target  ranges  eloser  than  1  km. 


Figure  44  displays  plots  of  the  power-in- the-bueket  data  for  the  same  study  as 
Figure  43.  Similar  to  previous  preeipitation  results,  power- in-the-bucket  inereases 
proportionally  with  target  proximity  to  the  source.  However,  a  100  kW  HEL  only 
becomes  effective  as  an  offensive  and  defensive  weapon  in  such  heavy  rains  for  target 
proximity  of  less  than  1000  m.  Otherwise,  the  actual  power  delivered  to  the  target  is 
much  too  low  to  inflict  significant  damage,  although  soft  kills  or  blinding  may  still  be 
achievable  at  longer  ranges. 


Standard  Model  •  Heavy  Rain  (25.0  mm/hr)  Target  Range  Study  •  PIB 

3000  m  2000  m  1000  m 

10®  1 - - - 


1  2  3  4  5 

Wavelength  (ijm) 


Figure  44.  Standard  model  heavy  rain  target  range  study:  Power-in- the-bueket 
results  increase  in  laser  performance  for  target  ranges  closer  than  1  km. 
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3,  Strong  Haze  Target  Range  Study 

Figure  45  displays  plots  of  the  laser  irradiance  results  for  the  standard  model  in 
the  presence  of  strong  haze  (5  km  meteorological  range),  within  which  the  proximity  to 
the  target  has  been  varied  from  1000  m  (indicated  by  dark  blue)  to  5000  m  (indicated  by 
light  blue).  In  the  presence  of  strong  haze  a  target  engagement  is  likely  to  be  successful  at 
a  range  of  5  km  or  less,  as  the  melt  threshold  is  exceeded  at  those  ranges.  The  increase  in 
on-target  laser  irradiance  is  dramatic;  as  the  target  approaches  to  within  3000  m,  the  laser 
intensity  increases  by  more  than  an  order  of  magnitude. 


Standard  Model  -  Strong  Haze  (5  km  VIS)  Target  Range  Study 
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Figure  45.  Standard  model  strong  haze  target  range  study:  Irradiance  results 
indicate  successful  engagement  at  target  ranges  closer  than  5  km. 


Figure  46  displays  plots  of  the  power-in- the-bucket  data  for  the  same  study  as 
Figure  45.  Closing  the  target  in  to  the  laser  source  results  in  more  power  delivered  to  the 
bucket.  As  the  beam  propagates  a  shorter  distance,  the  hazy  atmosphere  attenuates  the 
laser  power  less,  and  thus  the  on-target  laser  spot  becomes  more  lethal. 
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Figure  46.  Standard  model  strong  haze  target  range  study:  Power- in- the-bueket 
results  indicate  increase  in  power  delivered  to  the  target  for  ranges 

closer  than  5  km. 

4,  Strong  Turbulence  Target  Range  Study 

Figure  47  displays  plots  of  the  laser  irradiance  results  for  the  standard  model  in 
the  presence  of  strong  turbulence  (C^  =10”’‘*m~^^^),  within  which  the  proximity  to  target 

has  been  varied  from  1000  m  (indicated  by  dark  blue)  to  5000  m  (indicated  by  light  blue). 
A  successful  target  engagement  is  likely  for  target  ranges  closer  than  5000  m,  as  the  melt 
threshold  is  exceeded  at  those  ranges.  Since  the  laser  beam  has  to  propagate  a  shorter 
distance  to  the  target,  the  beam  is  less  affected  by  the  turbulent  environment  and  the  on- 
target  spot  size  thus  becomes  smaller,  while  the  laser  irradiance  becomes  greater. 


Figure  47.  Standard  model  strong  turbulence  target  range  study:  Irradiance  results 

indicate  likelihood  for  successful  engagement  at  target  ranges  closer  than  5  km. 
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Figure  48  displays  plots  of  the  power-in-the-bueket  data  for  the  same  study  as 
Figure  47.  It  appears  that  as  the  target  closes  in  on  the  position  of  the  laser  source,  the 
power-in-the-bucket  increases  with  the  target  proximity  though  less  so  at  ranges  closer 
than  2000  m.  For  these  closer  ranges,  the  effects  of  turbulence  are  so  small  that  the 
power-in-the-bucket  cannot  increase  with  all  of  the  available  laser  power  already 
captured  within  it. 


Figure  48.  Standard  model  strong  turbulence  target  range  study:  Power-in-the- 
bucket  results  indicate  increase  in  power  delivered  to  target  for  ranges 

closer  than  5  km. 


D,  DWELL  TIME  STUDIES 

In  addition  to  laser  irradiance  and  power-in-the-bucket  studies,  the  dwell  time 
performance  metric  (discussed  in  Chapter  IV,  Section  D.2)  was  also  incorporated  within 
the  standard  model  to  give  additional  insight  on  laser  performance  for  various  weather 
conditions  and  target  ranges. 

1.  Example  Dwell  Time  Calculation 

Given  a  target  with  parameters  as  described  in  Section  B,  the  required  amount  of 
time  the  laser  spot  needs  to  remain  on  target  in  order  to  inflict  significant  damage  can  be 
calculated  by  following  the  general  procedure  as  discussed  in  Chapter  IV,  Section  D.2. 
From  Equation  (67),  the  conductive  power  loss  is  calculated  as 
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(78) 


L 

and  the  radiative  power  loss  from  Equation  (68)  is  given  as 

P,=cjsA[t:-T:)  =  11^.  (79) 

The  total  power  loss  to  the  surroundings  is  thus 

=2.4  kW.  (80) 

In  this  example,  the  radiative  power  loss  to  the  surroundings  is  rather  negligible. 

Assuming  that  the  laser  delivers  50  kW  of  power  to  the  target,  the  available 
melting  power  is  ealculated  as  P^^,,  =(50kW-2.4kW)»48  kW.  Further  assuming  that 
the  target  is  sufficiently  damaged  once  0.1  L  of  its  skin  has  been  melted,  the  energy 
required  to  melt  this  volume  «153  kJ.  Finally,  the  dwell  time  required  to  melt  the 

0.1  liters  of  aircraft  skin  can  be  estimated  by  Equation  (72)  as 

^n=-^~3sec.  (81) 

melt 

Of  course,  the  actual  dwell  times  depend  upon  the  specifics  of  the  target.  The 
following  studies  aim  to  determine  the  dwell  times  required  for  a  laser  with  the  attributes 
listed  in  Table  7,  except  with  a  fixed  1.064  pm  wavelength,  to  effectively  melt  0.1  liters 
of  one  millimeter-thick  aluminum.  The  specified  volume  implies  the  assumption  that  the 
beam  spot  will  wander  while  on-target  in  order  to  melt  the  requisite  volume  of  material. 
The  dwell  times  are  determined  for  various  haze,  precipitation  and  turbulence  conditions 
and  at  different  target  ranges  in  order  to  fully  quantify  laser  lethality. 

2.  Dwell  Time  Study  for  Variable  Visibility  Conditions 

Figure  49  displays  plots  of  the  estimated  dwell  times  required  for  this  1 .064  pm 
wavelength  laser  to  inflict  damage  to  a  target  in  various  meteorological  range  conditions. 
As  evidenced  by  the  plot,  on  clear  days  with  visibility  of  23  km,  the  required  engagement 
dwell  times  fall  below  five  seconds  with  target  ranges  of  5  km  or  less.  On  hazy  days  with 
visibility  of  five  km,  the  dwell  times  fall  below  five  seconds  for  target  ranges  of  one 
thousand  meters  or  closer.  Note  that  the  5  km  visibility  curve  diverges  strongly  at  a  range 
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of  about  three  kilometers  while  the  eurves  eorresponding  to  less  hazy  eonditions  do  so  at 
longer  ranges.  This  is  the  result  of  the  faet  that  at  longer  ranges,  the  eonduetive  and,  to  a 
lesser  extent  the  radiative,  losses  prevent  the  aluminum  from  reaehing  its  melting  point 
regardless  of  the  dwell  time. 


Laser  Dwell  Time  Study  -  Variable  Haze  (Meteorological  Range)  Conditions 


Figure  49.  Laser  dwell  times  for  varying  meteorologieal  range  eonditions; 

Results  indieate  reasonable  engagement  times  for  target  ranges  less  than 
5  km  (dependent  on  visibility). 

3,  Dwell  Time  Study  for  Variable  Precipitation  Conditions 

Figure  50  displays  plots  of  the  estimated  dwell  times  required  for  a  laser  to  infliet 
damage  to  a  target  in  various  preeipitation  eonditions.  Even  small  amounts  of 
preeipitation  greatly  affeet  dwell  times;  however,  target  engagements  within  a  few 
kilometers  range  remain  viable.  For  heavy  preeipitation,  engagements  are  only  viable  for 
a  target  range  of  less  than  one  kilometer. 
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Laser  Dwell  Time  Study  -  Variable  Precipitation  Conditions 


Figure  50.  Laser  dwell  times  for  varying  precipitation  conditions:  Results 
indicate  viable  engagement  times  in  rainy  conditions  at  target  ranges 

less  than  1  km. 

4,  Dwell  Time  Study  for  Variable  Turbulence  Conditions 

Figure  5 1  displays  plots  of  the  estimated  dwell  times  required  for  a  laser  to  inflict 
damage  to  a  target  in  various  turbulence  conditions.  Dwell  times  remain  reasonable  even 
in  a  very  turbulent  atmosphere  for  target  ranges  of  five  kilometers  or  less. 


Figure  5 1 .  Laser  dwell  times  for  varying  turbulence  conditions:  Results  indicate 
reasonable  target  engagement  times  for  ranges  less  than  5  km,  even  in 
the  presence  of  turbulence. 
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VII.  CONCLUSIONS 


A.  CONCLUSIONS  FOR  5  KM  TARGET  RANGE 

According  to  the  results  for  the  standard  model  studies  at  a  5  km  target  range,  the 
100  kW-class  HEL  is  most  likely  to  successfully  engage  slow-moving  targets  (such  as 
light  aircraft,  UAVs,  etc.)  only  in  mild  atmospheric  conditions.  Specifically,  this 
requirement  is  met  if  atmospherie  turbulenee  is  weak  to  moderate  (C^  =10“''’m“^^^  to 

lO^^^m^^^^),  the  meteorological  range  is  greater  than  10  km  visibility,  and  no  appreciable 
preeipitation  is  present. 

Furthermore,  it  appears  that  large  beam  director  sizes  increase  on-target  laser 
irradianee  in  only  weak  turbulenee  eonditions.  It  is  quite  possible  that  moderate  to  strong 
turbulence  would  be  present  on  some  days  while  in  a  maritime  environment.  So,  laser 
weapon  systems  installed  on  naval  vessels  or  aircraft  could  be  designed  in  a  more 
compact  manner  with  smaller  beam  directors,  since  larger  beam  directors  provide  little 
added  benefit  in  sueh  conditions. 

B.  CONCLUSIONS  FOR  TARGET  RANGES  CLOSER  THAN  5  KM 

According  to  the  results  for  the  standard  model  studies  at  variable  target  ranges, 
the  100  kW-class  HEL  can  be  effeetively  utilized  in  an  offensive  and  defensive  capability 
against  slow-moving  targets  (such  as  light  aircraft,  UAVs,  etc.)  in  a  greater  variety  of 
atmospheric  conditions  as  the  target  range  decreases  below  five  kilometers.  Since  the 
various  eontributions  to  on- target  laser  irradianee  (diffraetion,  jitter,  turbulence, 
extinetion  and  thermal  blooming)  are  strongly  dependent  on  target  range,  these 
eontributions  eonspire  to  greatly  inerease  laser  irradianee  at  shorter  distanees. 

Furthermore,  these  effects  also  aid  in  the  total  power  delivered  to  the  desired 
target  area,  the  bucket.  The  closer  a  target  is  to  the  laser  source,  the  lesser  impact  that 
atmospheric  attenuation  and  beam  dispersion  have  on  the  laser  beam  as  it  propagates  to 
the  objective.  From  these  results,  100  kW-class  HELs  could  be  effectively  employed 
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against  targets  in  hazy,  turbulent,  and  rainy  conditions  as  long  as  the  target  is  no  more 
than  three  kilometers  from  the  laser  source,  depending  on  the  severity  of  the  weather. 


C.  CONCLUSIONS  FOR  VARYING  DWELL  TIME  STUDIES 

According  to  the  dwell  time  studies  in  clear  weather,  the  100  kW-class  HEL  can 
successfully  engage  and  damage  a  slow-moving  target  from  a  range  of  five  kilometers  or 
closer  within  a  five-second  dwell  time.  For  target  ranges  from  5  km  to  10  km,  longer 
laser  dwell  times  (~15  seconds)  are  required  to  inflict  a  similar  amount  of  damage.  For 
target  engagements  at  longer  ranges  and  in  poor  weather  conditions,  required  laser  dwell 
times  rise  to  greater  than  10  seconds,  and  in  some  cases  become  unattainable  (in  the  case 
of  heavy  rain,  for  example).  However,  for  hazy,  turbulent  and  rainy  days,  laser  dwell 
times  are  still  reasonable  at  shorter  ranges  (<  3  km). 
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